{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "58ecd0f8-90fc-429c-bb36-1d5002ed3972",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "import os\n",
    "import matplotlib as mpl\n",
    "from scipy.stats import lognorm\n",
    "from scipy.stats import beta \n",
    "import scipy as scp\n",
    "import matplotlib as mpl\n",
    "\n",
    "mpl.rc('font',family = 'Times New Roman')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "edfad454-d769-4777-a452-2946e2ed7c7c",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "d = np.zeros(30)\n",
    "\n",
    "for i in range(len(d)):\n",
    "    disp = np.loadtxt('./Results/dispI_GM'+str(i+1)+'.out')\n",
    "    d[i] = max(abs(disp[:,1]))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "3c79e4d7-9bfa-47c3-9ffb-320086573219",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "d = np.sort(d)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "29206d0a-973b-4b29-bddd-0c993b20a1e9",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.0230545"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "np.median(d)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "0e3c9007-42c8-41d6-bfd3-94d0c4fedfca",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "pc = np.linspace(0,1,len(d))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "a872d1c9-2b2b-4007-bcf5-a8eef699b595",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([0.        , 0.03448276, 0.06896552, 0.10344828, 0.13793103,\n",
       "       0.17241379, 0.20689655, 0.24137931, 0.27586207, 0.31034483,\n",
       "       0.34482759, 0.37931034, 0.4137931 , 0.44827586, 0.48275862,\n",
       "       0.51724138, 0.55172414, 0.5862069 , 0.62068966, 0.65517241,\n",
       "       0.68965517, 0.72413793, 0.75862069, 0.79310345, 0.82758621,\n",
       "       0.86206897, 0.89655172, 0.93103448, 0.96551724, 1.        ])"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pc"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "e5a76820-8495-44cb-9834-52679a8b2500",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x2e5c31b7ed0>]"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAh8AAAGdCAYAAACyzRGfAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABAF0lEQVR4nO3deXhU5d3/8fdMyMKSDCQsSSACZVNEUQSUugUFDSqoCFosgktRllq1Wh9CtZSKCj+ptVqFal1BQcANtCpuoJYlKuICCoIokARZgjNJIAmZOb8/jgkMmSQzySwnM5/Xdc3VmZP7nNznGR/z8Zz7fL82wzAMRERERMLEHukJiIiISGxR+BAREZGwUvgQERGRsFL4EBERkbBS+BAREZGwUvgQERGRsFL4EBERkbBS+BAREZGwahbpCRzL4/FQUFBAcnIyNpst0tMRERERPxiGQXFxMZmZmdjtdV/bsFz4KCgoICsrK9LTEBERkQbYuXMnnTp1qnOM5cJHcnIyYE4+JSUlwrMRERERf7hcLrKysqr/jtfFcuGj6lZLSkqKwoeIiEgT48+SCS04FRERkbBS+BAREZGwUvgQERGRsFL4EBERkbBS+BAREZGwUvgQERGRsFL4EBERkbBS+BAREZGwslyRMREREQkNt8cgb3sRe4rLaJ+cxMCuqcTZw99HrUHh480332T69OksXryYLl26+ByzbNky3nvvPcrLyxk1ahRDhgxpzDxFRESkEd76upAZyzdR6Cyr3pbhSGL68N7k9MkI61wCDh979uyhsrKSTz75pNYx3377LTNnzmTdunUYhkH//v1Zvnw5HTt2bNRkRUREJHBvfV3IpAXrMY7ZvttZxqQF65k7tl9YA0jAaz7at2/PxRdfXOeYhx56iJycHGw2G3a7nUGDBjF37twGT1JEREQaxu0xmLF8U43gAVRvm7F8E26PrxGh0aAFp3Z73bu9//77dO7cufpzjx49WLVqlc+x5eXluFwur5eIiIgER972Iq9bLccygEJnGXnbi8I2p5A87ZKfn09qamr15+TkZAoKCnyOvf/++3E4HNWvrKysUExJREQkJu0prj14NGRcMIQkfNhsNpKSkqo/V1RUEB8f73Nsbm4uTqez+rVz585QTElERCQmtU9Oqn9QAOOCISSP2mZmZuJ0Oqs/FxcXk5mZ6XNsYmIiiYmJoZiGiIhIzBvYNZUMRxK7nWU+133YgHSH+dhtuITkysf555/Pli1bqj9v3bqVwYMHh+JXiYiISB3i7DamD+8NmEHjaFWfpw/vHdZ6Hw0KH4ZheP0vwJw5c9i8eTMAEydO5N133wWgsrKSvLw8JkyY0Ni5ioiISAPk9Mlg7th+pDu8b62kO5LC/pgtNOC2S0lJCfPnzwfg2Wef5fe//z1t27Zl4cKFdOnShV69etG3b1+uu+467rjjDioqKvjHP/5Benp60CcvIiIi/snpk8HQ3umWqHBqM46+fGEBLpcLh8OB0+kkJSUl0tMRERERPwTy91u9XURERH5hld4n0U7hQ0REBGv1Pol2IXnaRUREpCmp6n1ybCXQqt4nb31dGKGZRSeFDxERiWlW7H0S7RQ+REQkplmx90nIVJbDnm/NVwRpzYeIiMQ0K/Y+aRTDgOLdsP872Pcd7N/6y/9+Bz/vAMMDJwyHqxZEbIoKHyIiEtOs2PvELxWlRwWLowLG/m1QUVL7fgnJEJcQvnn6oPAhIiIxzYq9T6p53ODcCfu2HnUl45eA4cqvfT9bHLTpDGk9oG0PSOtuvtr2gFYdwBbZx4cVPkREJKZV9T6ZtGA9NvAKIGHrfXLoQM2AsW8rFH0P7vLa92uR9kvA6O4dNNp0hWaRvbpRF4UPERGJeVW9T46t85EezDoflRVw4IeaAWP/Vji4r/b94hIgtVvNgJHWHVpE4GpMECh8iIiIEKTeJ4YBJXt8L/Y88CMY7tr3Tc70HTBaHwf2uMafoIUofIiIiPwizm5jULe0+gdWHISibb4Xe5a7at8vvqXvgJHWHRJbBe9ELE7hQ0REwqLJ9U3xeMC1q2bA2LfV3F4bm928WuFrsWdyRsQXe1qBwoeIiISc5fumGIa5uHPHGvhxDRRuMANHZR21PZq38R0wUn8FzRLDNvWmSOFDRERCqqpvyrGPsVb1TZk7tl/4A4jHDbu/gh1rYcdq839Lfqo5zh5vholjA0ZaD2jpx+0Z8UnhQ0REQqa+vik2zL4pQ3unh/YWzOFDkP+ZeVVjxxrYmQcVxd5j4hIgsx90HgSdBkK7XtC6M8TpT2Ww6f+iIiISMoH0TfFroae/DhaZAWPHajNwFHwOnsPeYxJTIOt0OO4M6PxrM3jEW6yKaZRS+BARkZAJW98U564jVzV2rIE9m2qOaZVuXtU47tdm4OhwYtQ9wtpUKHyIiEjIhKRviscD+7YcuaqxYy04d9Qcl9bjyFWN4wZBmy560sQiFD5ERCRkgtI3pbICCr84sjB0xxqzHLnXgeIg4+QjVzWOGwSt2gXzVCSIFD5ERCRkGtQ3pbwEduUduY2y61OoPOR94GbNoVP/I1c1Og2IqSJdTZ3Ch4iIhFS9fVO6xMGm18yrGj+uNh+BPbYMefNUM2RU3UbJ6Atx8WE+EwkWm2EYvq6ERYzL5cLhcOB0OklJSYn0dEREJEjcHoO87/dT8tN3dC39kl8d/BL7zrVmMa9jOY77ZXHoL6+2PcFuD/+kxW+B/P3WlQ8RkShmiZLmhgH5nxG38RUGbVrmY3GoDdr3PipsnAGOTuGdo4SVwoeISJSKaEnzXwIHG18xb6k4dx75mT0eOvYzg0bnX0PWQLNUucQMhQ8RkSgUkZLmXoHjmCsc8S2hVw6ceDl0Ox8SWgT3d0uTovAhIhJlwlrS3DAgfz1sfLn2wNH7MugxFOKbN+53SdRQ+BARiTIhL2leFTg2vQIbX1PgkIApfIiIRJmQlDSvL3D0vNC8pdJ9iG6pSL0UPkREokzQSpobBhSsP7Jo9GdfgeMy6D5UgUMCovAhIhJlGlXSvDpwvAqbXj0mcLSAnjkKHNJoCh8iIlEm4JLmhmG2nN/4Si2Bo+qWigKHBIfCh4hIFKq3pPmJ6b+s4XjVvMrx849Hdq4KHL0vgx4XKHBI0Cl8iIhEqZw+GQztnX6kwmmrRAYm7SDum7nw7qsKHBIxCh8iIlEszm5jUMp++P55WPlKzcDR4wJzDUePCyChZcTmKbFF4UNExIIa3ZPF44Ytb0He4/D9yiPbmzU/8pSKAodEiMKHiIjFNKonS+k+WP8cfPrUUf1UbOZTKn2vUuAQS1D4EBGxkAb3ZMn/DPKegK9fBne5ua15KvQbB/2vhzadQz53EX8pfIiIWETAPVkOl5mPx+Y9btbmqJJxCgy8EfqMVHlzsSSFDxERi/C3J8uGr77itL0vm7dXDu43fxiXACeOhIEToONpYGtkwziREFL4EBGxiLp7rRicaf+a8XEr6PfK54DH3JzSCQZcD6eOg1btwjFNkUZT+BARsQhfvVZacZAr4j7imrh36G4vOPKDrueat1Z65kCc/lUuTYv+iRURsYije7J0YD83NnuDK+NW0spmXhEpMZJ4M24wI2/6K3Edjo/sZEUaQeFDRMQi4uw2Zg1OpuD1R7gi7kMSbG4AvvN05Dn3BbziPos5Y88krkM9j9uKWJzCh4iIFezdDB89yLlfLYFmZuhY6zmBRysv5SPPSWQ4mjPHnzofIk2AwoeISCQVfgEf/R02LaO6/2z3IbjPuh3D3YtRxWVMbkiFUxELU/gQEQmhWsuk78yDDx+A71YcGXz8JXDOHZB5KnHAoIjNWiS0FD5EREKkZpl0g0uStzKj9X9J27vO3GSzQ58r4Kw/QofeEZurSDgpfIiIhIB3mXSDwfYN3NzsFfod3gp7wWNrhv2UMXDWbZDWLcKzFQkvhQ8RkSA7Uibd4Bz7l9zRbDEn27cDUG7Es8g9mFdaXMFLw6/SOg6JSQofIiJBlre9iEzXFzyU8CKn278FoNRIZL57KE9WXsReWoPLHDeoW1pkJysSAQofIiLBtPsruq74My8lrgLMKx3PuYcyt3IERaR4Da27nLpI9FL4EBEJBlcBvDMdvlpMOlBp2FnsPpdHKkdSiO+rG77KqYvEAoUPEZHGOFwGa/4FHz0Ih0sB8Jx4BVdvGcwnxalVlTu82IB0h/nYrUgsUvgQEWkIw4DNb8LbuXDgB3Nb1ukwbDb2zFO5/utCPlmwHht4BZCq5aXTh/fWYlOJWfZIT0BEpMnZuxkWjIRFY8zgkZwBI5+A69+GzFMByOmTwdyx/Uh3eN9aSXckMXdsP5VJl5imKx8iIv4qc8LK2ZD3b/BUQlwCDPo9nH07JLaqMTynTwZDe6f7rnAqEsMCDh+lpaXceeedtGnThpKSEmbPnk1iYmKNcffddx8ejwe73c7Bgwe55557sNn0/3Ai0gR5PLDheXhvBpTuNbf1HAYX3ltvgbA4u02P04ocI+DbLpMmTWLIkCHMnDmTfv36kZubW2PM66+/zvbt27nrrruYNm0aBQUFLF68OCgTFhEJq5158J/zYNnvzeCR1gN++xJcvcgreLg9Bmu27ee1Dfms2bYft8fXUlMRgQCvfBQUFLBkyRIef/xxAIYNG8bEiROZMWMGycnJ1eO++eYbiouLqz83b94cp9MZpCmLiIRB8W5496/wxULzc0IyZP8fDLwJmiV4Da3ZwwUyHElMH95baztEfAjoysfKlStp27YtSUnmAqp27dqRkJBAXl6e17gRI0awbNkyli5dyoEDB9i3bx/XXHONz2OWl5fjcrm8XiIiEVNZDh8/BI+cdiR4nDIWbv4Mfn2zz+AxacF6r+ABsNtZxqQF63nr68IwTVyk6Qjoykd+fj6pqd7PpScnJ1NQUOC1rVevXrzwwgtcffXVnHvuubz00ks0b97c5zHvv/9+ZsyYEeC0RURCYMvb8NZUKPre/NyxPwz7f9DpNJ/Dj/RwqcnAfKx2xvJNDO2drkWmIkcJ6MqHzWarvupRpaKigvj4+Bpjy8rKWLZsGdu2bWPMmDFUVlb6PGZubi5Op7P6tXPnzkCmJCLSePu2wvOj4YUrzeDRsj1cNhdueKfW4AFmb5Zjr3gczQAKnWXkbS8KwaRFmq6AwkdmZmaNtRslJSVkZmZ6bVu5ciWffvopQ4YMYeXKlXz++ec8+OCDPo+ZmJhISkqK10tEJCzKXLDibnjsDPhuBdjj4dd/MG+xnHI12Ov+V6S/vVnUw0XEW0DhIzs7m127dlFRUQFQfbtl4MCBXuOWLl1Kjx49ADOwzJw5kw8//DAY8xURCY7Nb8Kjp8Pqh8FzGLoPhclr4YJ7IMm//wjytzeLeriIeAv4ykdOTg6rVpndGlesWMHkyZNJSkpizpw5bN68GYBTTjmFzz//vHo/m81WI6CIiEREyV5Yej0s/A0UF0CbLnD1Yhi7FNp2D+hQA7umkuFIorbVHDbMp17Uw0XEW8BFxubNm8fUqVNZt24dRUVFzJo1C4CFCxfSpUsXevXqxfXXX8+OHTuYNWsWbdu2ZefOnUybNi3okxcR8ZthwJeLzQWlh4rAZjerk2bnQkKLBh0yzm5j+vDeTFIPF5GA2AzDsFQlHJfLhcPhwOl0av2HiARH8W547few9R3zc4c+MOIR6NgvKIdXnQ+RwP5+q7eLiES3H1fDkmuh5CezF8u5d8KZt0Jczaf0Gko9XEQCo/AhItHJMGDtXFhxFxhuaN8bRj0N7Y8Pya9TDxcR/yl8iEj0KS+BZTfDxpfNz31GwYiHIaFl9RC3x9CVCpEIUfgQkeiybyu8OBb2fgP2ZnDBvXD6TXBUV22t0RCJrIC72oqIWNY3y+HxbDN4tOoA41+HMybWCB7qxSISWQofItL0uSvNDrQvjoWKYjju13DTh9B5kPewenqxgNmLxe2x1EOAIlFH4UNEmrbSfbBgJHz8D/PzGVNg/DJITq8xVL1YRKxBaz5EpOna9RksHgeuXRDfEi59BPpcUetw9WIRsQaFDxFpegwDPnsG3rwT3BWQ1h2uWgDtT6hzN/ViEbEGhQ8RaVoOH4I37oANC8zPx18Cl831qxlcVS+W3c4yn+s+bEC6erGIhJzWfIhI03HgB3jyAjN42Oww5K/mFQ8/u9BW9WIBajSDUy8WkfBR+BCRpuG7d+Df58LuL6FFGlzzCpx1m9djtP7I6ZPB3LH9SHd431pJdyQxd2w/1fkQCQPddhERa/N44MMHYOX9gAEdT4MrnwNHpwYfUr1YRCJL4UNErOvQAXj5JvjubfPzadfBsNnQLDGgw9RWSl29WEQiQ+FDRKyp8EtYfI25zqNZElz8IJz624APo1LqItajNR8iYj0bFsKTQ83g0fo4uGFFg4OHSqmLWI/Ch4hYR2UFvHE7vDoRKsug+xC4cRVk9A34UCqlLmJdCh8iYg3OfHjmIvjkP+bnc6fC1UugRcNqbqiUuoh1ac2HiETe9g9hyXVwcB8kOWDkE9DzwkYdUqXURaxL4UNEIscwYPXDZkdawwMdToKr5kNq10YfWqXURaxL4UNEIqPMBa9Nhm+Wm5/7jjGfaEloEZTDq5S6iHVpzYeIhN/ezfCf883gYY+Hi/9u9mcJUvAAlVIXsTKFDxEJr42vwBPnwb4tkJwJ170JA34XcJl0f6iUuog16baLiISHuxLenQ5r/mV+7nI2jHoaWrUL6a9VKXUR61H4EJHQK9ljPs3y48fm5zNvgfP+AnHh+VeQSqmLWIvCh4gErLZeKT7tWAdLxkNxISS0gsseg96XhnfCImIpCh8iEhC/e6UYBuQ9AW9PA89haNsLrloA7XpGYNYiYiVacCoifvO7V0rFQXjlJnjzT2bw6H0ZTHhfwUNEAF35EBE/1dcrxYbZK2Vo+kHiloyDn74GWxwM/RsMmhKSp1lEpGlS+BARv/jTK6V38f8w/j0eDhdDy3Yw+hnoclbY5igiTYPCh4j4pa4eKHY83NpsKX9o9iocBrJON4NHSma4piciTYjCh4j4pbYeKK0p5uH4f3FO3FcAFPYaT8boOdAsIZzTE5EmRAtORcQvVb1Sjl650cf2Pa8n/plz4r7ikJHA9Ga30P6qfyp4iEidFD5ExC/H9koZHbeSlxJm0Mm2jx88HRhZ8TcGXTZJlUNFpF667SIifsvpk8Hc357KT6/+mfHulwF4x92P/9f8Nm6/coB6pYiIXxQ+RMR/Hjc5P8yGX4LHt8dPoVX/23jrV211xUNE/KbwISL+qawwC4dtfBnDZueLvtM5dNI1atImIgFT+BCR+lUchMXjYOs7HKYZfyifwptre8Datb5Lq4uI1EELTkWkbmVOWDAStr7DISOBGypu503P6dU/rlFaXUSkHgofIlK7kr3wzCWwYw3FtGBsRS4fevp6Dakqtz5j+SbcHl/F10VEvCl8iIhvzl3wdA7s/pKKpDSuKr+Lz4xePocaQKGzjLztReGdo4g0SQofIlLTzzvh6Ytg/1ZwZPHhWfPZZHSpd7e6SrCLiFRR+BARbz/vhGcuhp9/hDZd4fq3aJlxvF+71laCXUTkaAofInKEcxc8e8mR4HHtG+Do5LO0+tFsQIYjiYFdU8M5WxFpohQ+RMTk3GVe8TjwA7TpAte+Do6OQM3S6ker+jx9eG/V+xARvyh8iMgvweOSo4KHecXjaDl9Mpg7th/pDu9bK+mOJOaO7ac6HyLiNxUZE4l1zvxfgsf2WoNHlZw+GQztnU7e9iL2FJfRPjlJFU5FJGAKHyKxzJn/y62W7dC6M4x/vdbgUSXObmNQt7QwTVBEopHCh0ischWYi0urgse1b0DrLADcHkNXN0QkZBQ+RGKRq8C84lH0/S/B4/Xq4PHW14XMWL6JQueRmh3q3yIiwaQFpyKxxlVgrvEo+h5aH/dL8DgOMIPHpAXrvYIHqH+LiASXwodILKkOHtt+CR5vVAcPt8dgxvJN+OrOov4tIhJMCh8iscJVeCR4OI4zF5f+EjwA8rYX1bjicTT1bxGRYFH4EIkFrkJzcWlV8Lj2dWjT2WuIv31Z1L9FRBpL4UMk2hXvNoPH/q21Bg/wvy+L+reISGMpfIhEs+Ld5lMtv3SnrS14AOrfIiJhE/CjtqWlpdx55520adOGkpISZs+eTWJios+xBw8e5N///jdpaWn07NmTM844o9ETFhE/Fe8213j4ETzgSP+WSQvWYwOvhafq3yIiwRTwlY9JkyYxZMgQZs6cSb9+/cjNzfU5rqioiNGjR3PZZZcxbtw4BQ+RcCr+CZ4dDvu/g5ROMH65WTq9HurfIiLhYDMMw+/n5goKCujWrRsHDhwgKSmJvXv30rlzZ3766SeSk5O9xl5wwQVMmzaN7OzsgCbkcrlwOBw4nU5SUlIC2ldEgINF5q2WPZvM4HHt65DaNaBDqMKpiAQqkL/fAV35WLlyJW3btiUpyfyvonbt2pGQkEBeXp7XuNdff53vvvuOvLw8LrroInJzczl8+LDPY5aXl+NyubxeItJAZS5YcIUZPFqlw7XLAw4ecKR/y6WndGRQtzQFDxEJqoDCR35+Pqmp3ovNkpOTKSgo8Nr2/PPPc/rpp3Prrbcyf/58nn/+ee6++26fx7z//vtxOBzVr6ysrABPQUQAOFwGi66GgvUcTmzD+wMfZ80Bh4qCiYjlBBQ+bDZb9VWPKhUVFcTHx3tt27hxI2eddRYJCQmkpaUxYcIEnnvuOZ/HzM3Nxel0Vr927twZ4CmICB43vDwBfviIUpIY6bqd6/9bwpgn1nLW7PdVFl1ELCWg8JGZmYnT6fTaVlJSQmZmpte2yspK3G539eeTTz6ZoiLfVRETExNJSUnxeolIAAwD/nsHfLOMcqMZEyr+yFfGr6p/rL4sImI1AYWP7Oxsdu3aRUVFBUD17ZaBAwd6jTv55JP57rvvqj83a9aM448/vrFzFRFfVv0/+PQpPNj44+HJrPb08fqx+rKIiNUEfOUjJyeHVatWAbBixQomT55MUlISc+bMYfPmzQDccsstvPHGG5SXlwOwevVqbr755iBPXUT49ClYeR8Afz08jjc8vh9pV18WEbGSgIuMzZs3j6lTp7Ju3TqKioqYNWsWAAsXLqRLly706tWLQYMGcc8993DzzTfTq1cvkpKSuOGGG4I+eZGYtmkZvHE7AJt7TuS5L8+pdxf1ZRERKwiozkc4qM6HiB+2fwQLRoK7AvqNZ03vvzDmP+vq3W3hhDMY1C0tDBMUkVgTyN/vgK98iEiE7f7KfKTWXQHHXwIXP8hAWxwZjiR2O8vw9V8TNswqperLIiJWoMZyIk1J0XaziFi5C477NVzxH4hrVt2XBajRGE59WUTEahQ+RJqKkr3mrZaSn6D9iTBmIcQ3r/6x+rKISFOh2y4iTUF5MTw/Coq+h9bHwdiXoHnrGsNy+mQwtHe6+rKIiKUpfIhYXWUFvDgWCjdAizQY+wqk1LyKcWwzuEtOzlToEBFLUvgQsTKPB16dCN+vhPiW8Nsl0LZ7jWFvfV3IjOWbKHQeeZQ2w5HE9OG9dbtFRCxHaz5ErMow4O1c+PolsMfDVfOh42k1hr31dSGTFqz3Ch6gsuoiYl0KHyJWtfoRWDfPfH/5POh+fo0hbo/BjOWbfD5eq7LqImJVCh8iVvT1y/DO3eb7C+6Fk0b5HJa3vajGFY+jqay6iFiRwoeI1fy4Bl6ZaL4/fSIMmlLrUH/LpausuohYicKHiJXs3waLxoC73KxeeuF9YKv9iZX2yUm1/qwh40REwkHhQ8QqSvebtTwOHTAXlo58Auxxde4ysGsqGY6kGlVNq9gwn3pRWXURsRKFDxErcFfCkvFHioiNWQQJLerdTWXVRaQpUvgQsYJ3p8MPH0FCK7h6MbRq7/euKqsuIk2NioyJRNpXS2HNv8z3l82F9icEfAiVVReRpkThQySSdn8Fr/3efH/WH6H3iAYfKs5uY1C3tCBNTEQkdBQ+RCLlYBEs+i1UHoJu58N5d3n9+NheLbqSISLRQuFDJBI8bnjpBvj5R2jTBa74j9eTLerVIiLRTAtORSLh/Xtg2/sQ3wKueh5aHHkUVr1aRCTaKXyIhNvGV+Hjf5jvRzwC6X2qf6ReLSISCxQ+RMJpzzfw6mTz/aDf1+jZol4tIhILFD5EwuXQz+YC08Ol0PUcGDKjxhD1ahGRWKDwIRIOHg+8fCMUbQNHFox6BuJqrvdWrxYRiQUKHyLhsGoWfPc2NEuCqxZAS9/1ONSrRURigcKHSKh9+19YNdt8f8lDkHlKrUPVq0VEYoHCh0go7fvOvN0CMPAmOGVMvbuoV4uIRDsVGRMJlTIXLLoaKorhuF/Dhff6vat6tYhINFP4EAkFjwdenQT7tkByJlz5LMTF17ubSqqLSCxQ+BAJhY8fhG9fh7gEuGo+tGpf7y4qqS4isUJrPkSC7bt34P2Z5vuL5kCn/vXuopLqIhJLFD5Egmn/NrNhHAacdh2cNr7eXVRSXURijcKHSLCUl8CLY6HMCZ0GwLDZfu2mkuoiEmsUPkSCwTBg2e9hzyZo1QGunA/NEv3aVSXVRSTWKHyIBMPqR2DjK2BvBqOfhRT/F4iqpLqIxBqFD5HG2vYBvDvdfJ8zCzoPCmh3lVQXkVij8CHSGAd+hKXXg+GBU34LA34X8CFUUl1EYo3Ch0hDVRyEF38Lh4og81S4+EGwNSwgqKS6iMQSFRkTaQjDgNdvhd1fQYu2Zqfa+MatyVBJdRGJFQofIg2x7t/w5Ytgi4PRz4CjU1AOG2e3MahbWlCOJSJiVQofIoH64WN4e5r5/oKZ0PXsoBxWfV1EJFYofIgEwrkLFo8Hww0njYYzJgXlsOrrIiKxRAtORfx1uAxevAYO7oMOJ8Hwhxu8wPRo6usiIrFG4UPEH4YB/70dCtZD8zbwmwWQ0KLRh1VfFxGJRQofIv749Cn4fAHY7DDqKWjTJSiHVV8XEYlFCh8i9dmxDt78P/P9+X+BbucF7dDq6yIisUjhQ6QurkJYfA14DkPvS+HMW4N6ePV1EZFYpPAhUpvKClgyHkp+gnYnwKWPBWWB6dHU10VEYpHCh0ht3poKO9dBogN+8zwktgr6r1BfFxGJRQofIr6snw+fPgnY4IonIK1byH6V+rqISKxRkTGRY+36DN74o/l+8DToeWHIf6X6uohILNGVD5Gjlew1F5i6K6DXxXD2HZGekYhI1NGVD5Eq7sOw5Fpw5UNaD7h8HtjDk89VXl1EYomufIhUWXE3/PgxJCTDb16ApJSw/FqVVxeRWKPwIQLwxYuwbq75/vJ50K5nWH6tyquLSCxS+BAp/AKW/8F8f86f4IRLwvarVV5dRGKRwofEttL9sGgsVJZB96GQnRvWX6/y6iISiwIOH6WlpUyZMoW77rqLW2+9lfLy8jrHL1q0iOzs7IbOTyR03JWw9Dpw7oA2Xc16Hva4sE5B5dVFJBYFHD4mTZrEkCFDmDlzJv369SM3t/b/UiwoKODee+9t1ARFQua9GbB9FcS3NBeYNm8T9imovLqIxKKAwkdBQQFLlixh2LBhAAwbNox58+ZRXFzsc/ysWbOYNGlS42cpEmxfvwyrHzbfX/ov6NA7ItNQeXURiUUBhY+VK1fStm1bkpLMS8Dt2rUjISGBvLy8GmOfeOIJxo4dS4sWLeo8Znl5OS6Xy+slElI/bYTXppjvz7wF+oyM6HRUXl1EYk1ARcby8/NJTfW+/JucnExBQYHXtq1bt+JyuRg4cCCbNm2q85j3338/M2bMCGQaIg136AAsuhoOH4RfZcN5f4n0jACVVxeR2BJQ+LDZbNVXPapUVFQQHx9f/dntdvPYY4/xwAMP+HXM3Nxc/vjHP1Z/drlcZGVlBTItEf943PDSBDjwA7Q+DkY9DXHWKfIbZ7cxqFtapKchIhJyAf2bNzMzE6fT6bWtpKSEzMzM6s+rV69m3rx5PPXUU4AZTioqKmjdujU///xzjWMmJiaSmJjYgKmL1M7tMWpeRfjgPtj6DjRLgqsWQAst4hQRiYSAwkd2djY33ngjFRUVJCQkVN9uGThwYPWYAQMGeN1qWbp0KUuXLmXRokVBmrJI3Xz1SflNqw3Mqpxjfhj+MGT0jdDsREQkoAWnmZmZ5OTksGrVKgBWrFjB5MmTSUpKYs6cOWzevJmkpCS6dOlS/apaoNqlS5dQzF/Ei68+Kd1s+dx12Hyy5Yfu46DvVZGanoiI0ICutvPmzWPq1KmsW7eOoqIiZs2aBcDChQvp0qULvXr1CvokRfzhq09KMgd5PP5BWtnKWOs5gTt2XMIqj6GFnCIiEWQzDMNSHatcLhcOhwOn00lKSni6ikp0WLNtP2OeWFv92YaHx+P/wdC4zygwUhlefi/7cbBwwhla2CkiEmSB/P1WbxeJGsf2P7k57lWGxn1GuRHPxIrb2I/D5zgREQkvhQ+JGkf3PznL/hW3NnsJgLsqr+NLo5vPcSIiEn4KHxI1qvqkdOAAD8U/it1m8ELlYJa4swH1SRERsQrrVFgSaaQ4u42/XtwTx9I/09bmYpOnMzMqxwPqkyIiYiUKHxJVLtz7DNi/oZQkphz+A+UkAGaflOnDe6tPioiIBSh8SPTY+i589HcAmo98lPtaDlafFBERC1L4kOjgKoCXbwQM6H899pNHMSjScxIREZ+04FSaPnclxtLr4eB+fnYcz9qef8LtsVT5GhEROYrChzR525ZMw7ZjDcVGcy7bM4HfPPU5Z81+n7e+Loz01ERExAeFD2nSPn13Md2+/TcAuYd/xw+GuaB0t7OMSQvWK4CIiFiQwoc0We6f8+n28e0AzK8cwuueI6s8qm66zFi+SbdgREQsRuFDmiZ3JaUvjKcNLjZ6OjOzcmyNIQZQ6Cwjb3tR+OcnIiK1UviQpumDe0nZ8wnFRnMmH76lup6HL+rlIiJiLQof0vR89y58/CAAUw9P4Ecjvc7h6uUiImItCh/StDjz4ZUbAfD0/x3rk7OprXSYermIiFiTwoc0He5KeOkGOLgf0k/GfuG9TB/eG6BGAFEvFxER61L4kKbjg3thxxpISIbRz0B8Ejl9Mpg7th/pDu9bK+mOJOaO7adeLiIiFqTy6tI0HLXOg0sfgbRu1T/K6ZPB0N7p5G0vUi8XEZEmQOFDLMntMarDRCf7Afq9OcG8lTJgApx4eY3xcXYbg7qlhX2eIiISOIUPsZy3vi5kxvJNFDrLiMPNooR7sNmLcLbujeOCmZGenoiINJLWfIilvPV1IZMWrKfQadbmuL3ZEgbYt+AymjPipwm8tflAhGcoIiKNpfAhluH2GMxYvqm6NHq2fQOTmy0DzHoeO4wOKpcuIhIFFD7EMvK2F1Vf8chkHw/GPwbAs5VD+a/nDJVLFxGJEgofYhlVZdDjqeTRhIdJtZXwpacr9x7Tt0Xl0kVEmjaFD7GMqjLouc1e4FT7VpxGCyYfvoUK4n2OExGRpklPu4hlDOyaym9bref6yrcA+OPhSewy2lf/3IZZPEzl0kVEmjZd+RDLiCvaxgzmAjC3cjjveU6r/pnKpYuIRA+FD7GGioOweBzNKkspatuf51uM8/qxyqWLiEQP3XYRa/jvn2DPRmjZntTxz7OqZQeVSxcRiVIKHxJ56+fDhgVgs8OoJyE5nThQuXQRkSil8CFBc3Q/Fr+vVuz+Cv57h/l+8DToek7oJyoiIhGl8CFBcXQ/lioZjiSmD+9d+zqNMicsHgeVZdB9KJx1e5hmKyIikaQFp9Jox/ZjqbLbWcakBet56+vCmjsZBrz2eyj6HhxZMPJxsOsfRxGRWKB/20ujHNuP5WhV23z2Y1k3D75ZBvZ4GP0MtFDtDhGRWKHwIY1ydD8WX3z2Y9mZByvuMt9feC906h/aSYqIiKUofEij+NtnpXpc6X5Yci14KuHEy2HgjaGbnIiIWJLChzSKv31W2icngccDL08AVz6kdYcRj4BNtTtERGKNwoc0ysCuqWQ4kqgtQtgwn3oZ2DUVPpoD296DZs3hyucgMTmcUxUREYtQ+JBGibPbmD68N0CNAOLVj2X7SvjgPnPDJQ9ChxPDNUUREbEYhQ9ptJw+Gcwd2490h/ctmOp+LMcZ8NLvAANOvQZOuToyExUREUtQkTEJipw+GQztnV6zwqlRCc9cAgf3QYeT4KIHIj1VERGJMIUPCVhtZdTj7Laa/VhWzICdayExBa58FuKbR2bSIiJiGQofEpCAyqh/8zqsfsR8f+mjkNYtjDMVERGr0poP8VtAZdSLtsOrk833Z0yB3iPCOFMREbEyhQ/xS0Bl1A+XmQ3jyp3QaSAMnRHOqYqIiMUpfIhfAiqj/tb/we4voUUajH4a4uLDN1EREbE8rfkQv/hbRj1h42JY/wxgg5FPgKNTSOclIiJNj658iF/8KaPey7aDU7745RbLuf8H3c8P8axERKQpUvgQv9RXRt1BCU8m/oM49yH41WA4986wzk9ERJoOhQ/xS11l1OPw8FD8o3TiJ2h9HIx6Cuxx4Z+kiIg0CQof4rfayqjf1fJVBsd9Ac2S4KoF0CI1QjMUEZGmQAtOJSDHllHv9fNHHL9yqfnD4Q9DRt/ITlBERCxP4UMCVl1Gfd938N87zI2nT4S+V0V2YiIi0iQofMSY2vqyBKy8GBb9FiqK4bhfwwUzgz9ZERGJSgofMSSgvix18Xjg1UmwbzMkZ8DoZ1RITERE/KYFpzEioL4s9Vk1C75ZDvZ4uHI+JHcI8mxFRCSaBRw+SktLmTJlCnfddRe33nor5eXlNcY4nU5Gjx5NSkoKp556KmvXrg3KZKVhAurLUp8vF8Oq2eb7S/4BWQOCNU0REYkRAYePSZMmMWTIEGbOnEm/fv3Izc2tMWb27NmMHDmSDz74gKysLC699FJKS0uDMmEJXEB9WeqyYx28NsV8f+Yt0O+a4E1SRERiRkDho6CggCVLljBs2DAAhg0bxrx58yguLvYaN2TIEMaMGcNpp53GggULOHDgAJs2bQrerCUg/vZlqXPcgR9h0dXgroBeF8P5fw3O5EREJOYEFD5WrlxJ27ZtSUoyi0y1a9eOhIQE8vLyvMadd9551e9TUlJISUmhUyffDcbKy8txuVxeLwkuf/qy1DmuzAULfwMH90H6STDycbBruZCIiDRMQH9B8vPzSU31rl6ZnJxMQUFBrfts2bKF7OxsMjJ8P01x//3343A4ql9ZWVmBTEn8UF9fFhvmUy8Du/qoTOquhKXXw55N0CodxrwIia1COV0REYlyAYUPm81WfdWjSkVFBfHxtT9m+eijj/LAAw/U+vPc3FycTmf1a+fOnYFMSfxQV1+Wqs/Th/f2Xe9jxV2w9R1o1hzGLARHx5DOVUREol9A4SMzMxOn0+m1raSkhMzMTJ/j3377bc477zy6du1a6zETExOrb81UvST4auvLku5IYu7Yfr7rfHzyJKyba76/fB507BeGmYqISLSzGYbhx/OVpoKCAnr27ElRUREJCQkUFBTQo0cP9u/fX+OKyMaNG1m9ejUTJkwIaEIulwuHw4HT6VQQCQG/K5xu+wAWXAGGG867G865I/yTFRGRJiOQv98BX/nIyclh1apVAKxYsYLJkyeTlJTEnDlz2Lx5MwDbtm3jiSeeYOjQofzwww988cUXzJ07t4GnI8FU1Zfl0lM6Mqhbmu/gsXcLLB5vBo+TfwNn3x7+iYqISNQKuLz6vHnzmDp1KuvWraOoqIhZs2YBsHDhQrp06UJKSgrZ2dns2rWLf/7zn9X7LVq0KHizllo1undL6X544Uood0LWGTDiYbA1oPeLiIhILQK67RIOuu3ScI3u3VJZDs9dBjtWQ+vOMOF9aNk2dBMWEZGoEbLbLmJdje7dYhjw+m1m8EhMgatfVPAQEZGQUPiIAkHp3fK/f8KG58Fmh9FPQ/sTQjFVERERhY9o0OjeLd8sh3f/ar7PmQ3dhwR9jiIiIlUUPqJAo3q3FGyAl28EDBgwAU6/MahzExEROZbCRxRocO8WV6HZs+XwQeh2HuTMCsHsREREvCl8RIEG9W6pOGgGj+JCaHc8jH4G4gJ+8lpERCRgCh9RIODeLR4PvHITFG6AFmkwZhEkOcI1XRERiXEKH1EioN4tH8yEb5ZBXAJc9Tyk1t57R0REJNh0nT2K5PTJYGjv9LornOY9AR/93Xw//GHoPCgykxURkZil8GFxgZZLr+rd4tOXS+C/fzLfnzsVThkTghmLiIjUTeHDwhpdLv1oW96GVycCBgy8EbKnBneyIiIiftKaD4tqdLn0o/24GhaPA08lnHSlWUhMzeJERCRCFD4sKCjl0qsUfgEvXAWVZdAzBy57DOz62kVEJHL0V8iCGl0uvcq+rTB/JJS7oPOZv9TyiA/qXEVERAKl8GFBjSqXXsWZD/Mvg4P7IP1kGLMQ4psHZ4IiIiKNoPBhQQ0ul16ldD/MvxycOyGtO4x9WUXERETEMhQ+LKhB5dKrlDnh+VGwbzOkdIRrXoFW7UI5XRERkYAofFhQwOXSqxwsgucuhYL10DzVDB6tjwv5fEVERAKh8GFRAZVLByjdB8+OgILPzeAx7jVo1yuMMxYREfGPioxZmF/l0gGKf4LnRsDeb6FlezN4dOgdmUmLiIjUQ+HD4uoslw7mUy3PjYD9WyE5A8Yvh7Y9wjdBERGRACl8WFi9fV0O/AjPDoeffwRHFoxfBqm/ityERURE/KDwYVH19nXZv81cXOrcCW26mFc8tLhURESaAIUPC6rq63Js8fSqvi7PjWjN2WtugOJCs47H+OWQkhmRuYqIiARK4cNi6uvr0su2kxNXTAKc0O4Ec3Fpcocwz1JERKTh9KitxdTV1+VE23YWJtxDKk5K2/SGa99Q8BARkSZH4cNiauvXcoptKy8k3EuqrYQNnl+xatB/oGUdT8GIiIhYlMKHxfjq19Lf9i3zE+7HYTvIp56eXFMxjTZp6RGYnYiISOMpfFjMsX1dLrB/woKE+0m2HWKNuzfjK6bSypHqu6+LiIhIE6DwYTFH+roY/C7uDebFP0SS7TDvuk/l+sN/4iBJvvu6iIiINBF62sWCck5oy4cnvErW9iUAPF95Pn+pvJb2jpZH6nyIiIg0UQofVnOwCJaMJ2v7hxjY+KH/n2nVcQwLUpr77usiIiLSxCh8WEBVGfXSwm85M28KzV3bIaEVtiuepGuvHLpGeoIiIiJBpPARYVVl1DsXr2de/D9obiulkLZsO+dJzuqVHenpiYiIBJ0WnEaQWUb9My4seZX58ffT2lbK557ujCi7h2teL+WtrwsjPUUREZGgU/iIELfH4P8t+4x/xv+Lv8Y/R7zNzWvuX/ObirvYiwOAGcs34fb4KrQuIiLSdOm2S4R8ueET5pXdSc+4fA4bcdxb+VuecV8Iv1T4MIBCZxl524sY1E2VTEVEJHoofISbcxd89Hf6fvYcdnslu402TKn4A58ZvXwOr63cuoiISFOl8BEurkL4+EH47BlwV2AHPnD35U+HJ7Lvl9ssvvgqty4iItKUKXyEWvFP8L+H4JMnwV1ubut8Fu7sXKYtKmd/LR1sbUC6I0ll1EVEJOoofIRKyd4joaPykLkt6ww478/Q9RzigOnDC5m0YD02zDUeVarKiKmMuoiIRCOFj2Ar3Q+rH4a8x+HwQXNbpwEweBr8ajDYjoSJnD4ZzB3bjxnLN1F41BWQdEeSyqiLiEjUUvgIloNFsOZRWDcPKkrMbZmnwuA/Q/chXqHjaDl9MhjaO5287UXsKS6jfXKSyqiLiEhUU/horDInrJ1rBo9yl7kt/SQzdPTMqTV0HC3ObtPjtCIiEjMUPhqqvNi8yrH6ETOAALQ/EQbnwvGX1Bk6qnq56EqHiIjEIoWPQJWXmOs5Vj8Mhw6Y29odD9lT4YRLwV530diqXi5Hr/HI0BoPERGJIQof/qo4CJ/8x3yC5eB+c1tad8jOhRMvB3tcvYcwe7ms59iC6budZUxasJ65Y/spgIiISNRT+KjP4UPw6dPw8T+gdI+5rU1X80pHn1EQ59//Cd0egxnLN9UIHmA+ZmvD7OUytHe6bsGIiEhUU/ioTWU5fPYsfPR3KNltbmt9HJz7f3Dyb/wOHVXythd53Wo5lnq5iIhIrFD4OFZlBXw+3wwdrnxzmyMLzrkD+l4NzRIadFh/e7Sol4uIiEQ7hY8q7sOw4QX4cA44d5jbkjPhnNvh1GugWWKjDu9vjxb1chERkWin8OGuhC9fhFWz4ecfzW2tOsDZt0O/8RAfnDAwsGsqGY4kdjvLfK77UC8XERGJFbEbPjxu+GoprJoFRd+b21q2g7Nug/7XQ3zzoP66OLuN6cN7q5eLiIjEvNgLHx4PbHzZvNKxb4u5rXkqnHUrDPgdJLQM2a9WLxcREZFYCh8eD3yzDFbOgr3fmNuSWsOZf4CBN0JiclimoV4uIiIS62ImfHjW/Rv721MBqIxPxn7mzdjPmARJKWGfi3q5iIhILAs4fJSWlnLnnXfSpk0bSkpKmD17NomJNZ8EWbZsGe+99x7l5eWMGjWKIUOGBGXCDfHW14U8+EEmTxtpLHWfy5Nlw2i5Jo3pbUvJ6RP+8CEiIhLLbIZh+Hr4olbjxo3j8ssv5/LLL+e5555jw4YNPPjgg15jvv32W8aNG8e6deswDIP+/fuzfPlyOnbsWO/xXS4XDocDp9NJSkrjg8HRJc3jcOPGLINedZNDJc1FREQaL5C/33V3QTtGQUEBS5YsYdiwYQAMGzaMefPmUVxc7DXuoYceIicnB5vNht1uZ9CgQcydOzfA02i8Y0uaVwUPOPK0yYzlm3B7AspfIiIi0ggBhY+VK1fStm1bkpLM2hft2rUjISGBvLw8r3Hvv/8+nTt3rv7co0cPVq1a5fOY5eXluFwur1ewBFLSXERERMIjoPCRn59Paqp3Eazk5GQKCgrqHOdrTJX7778fh8NR/crKygpkSnVSSXMRERHrCSh82Gy26qseVSoqKoiPj69znK8xVXJzc3E6ndWvnTt3BjKlOqmkuYiIiPUE9LRLZmYmTqfTa1tJSQmZmZl1jisuLq4xpkpiYqLPp2WCQSXNRURErCegKx/Z2dns2rWLiooKgOpbKQMHDvQad/7557Nly5bqz1u3bmXw4MGNnWvAqkqaw5GnW6qopLmIiEhkBBQ+MjMzycnJqV48umLFCiZPnkxSUhJz5sxh8+bNAEycOJF3330XgMrKSvLy8pgwYUKQp+6fqpLm6Q7vWyvpjiQ9ZisiIhIBAdf52LdvH1OnTqVLly4UFRUxa9YsEhISOO2008jNzWXUqFEAPP3002zcuJGKigouv/xyv698BLvORxW3x1BJcxERkRAJ5O93wOEj1EIVPkRERCR0QlZkTERERKSxFD5EREQkrBQ+REREJKwUPkRERCSsFD5EREQkrBQ+REREJKwUPkRERCSsFD5EREQkrBQ+REREJKwC6mobDlUFV10uV4RnIiIiIv6q+rvtT+F0y4WP4uJiALKysiI8ExEREQlUcXExDoejzjGW6+3i8XgoKCggOTkZm81ajd9cLhdZWVns3LkzJvrOxNL5xtK5QmydbyydK+h8o5nVz9UwDIqLi8nMzMRur3tVh+WufNjtdjp16hTpadQpJSXFkl98qMTS+cbSuUJsnW8snSvofKOZlc+1viseVbTgVERERMJK4UNERETCSuEjAImJiUyfPp3ExMRITyUsYul8Y+lcIbbON5bOFXS+0SyaztVyC05FREQkuunKh4iIiISVwoeIiIiElcKHiIiIhJXCh4iIiIRVzIaP0tJSpkyZwl133cWtt95KeXm5z3HLli3jlltuYeLEibz77rteP9u3bx933nknU6ZMCWi/cAv1uZaWlpKamorNZsNms/HKK6+E5Dz81djzdTqdjB49mpSUFE499VTWrl3r136REurztdL329hzdbvdTJw4kZSUFE466SQ+//xzv/aLlFCfbzR9t0dbtGgR2dnZAe8XTqE+Xyt9tz4ZMeqaa64xXn75ZcMwDOPZZ581brvtthpjvvnmG2PAgAGGx+Mx3G63ceqppxq7du2q/vlnn31mjBgxwhg/fnxA+4VbKM/VMAzj4YcfNhYvXmy88847xjvvvGNUVlaG7Fz80djzzc3NNV544QXj008/NYYPH260b9/eKCkpqXe/SAnl+RqGtb7fxp7rggULjDVr1hh79uwxLrroIuPMM8/0a79ICeX5GkZ0fbdV8vPzjT59+hjnnntuQPuFWyjP1zCs9d36EpPhIz8/30hKSjIOHTpkGIZh7Nmzx2jevLnhcrm8xt10003G3XffXf158uTJxp///GevMXfffXeNP8j+7BcuoT7XyspK44ILLjC+/fbb0JxAgIJxvu+99171dqfTacTHxxt5eXn17hcJoT5fK32/wTjXgwcPVm9fvny5cc455/i1XySE+nyj7butcvPNNxuPPvqo1x/jaPxuq/g6Xyt9t7WJydsuK1eupG3btiQlJQHQrl07EhISyMvL8xr3/vvv07lz5+rPPXr0YNWqVV5jfDXP8We/cAn1ua5YsYL//e9/HH/88Vx44YXs2bMnBGfhv2Cc73nnnVe9vaqHQlW/ISt9txD687XS9xuMc23evHn19k2bNnHffff5tV8khPp8o+27BXjiiScYO3YsLVq0CGi/cAv1+Vrpu61NTIaP/Px8UlNTvbYlJydTUFBQ5zhfY/w5vr/7hUKoz3XYsGEUFxezatUq8vPzGT58OB6PJziTb4Bgn++WLVvIzs4mIyMjoP3CJdTna6XvN1jnun//fv72t7/x6KOP4nQ6/d4v3EJ9vtH23W7duhWXy8XAgQPrPX40fLd1na+VvtvaxGT4sNls1YmzSkVFBfHx8XWO8zXGn+P7u18ohPpcq/Y955xz+OCDD9i6dWuNBYvhFOzzffTRR3nggQcC3i9cQn2+Vfta4fsN1rkmJydz0UUXMWDAAC6//HJ27tzp137hFurzrdo3Gr5bt9vNY489xq233urX8Zv6d1vf+Vbta4XvtjbNIj2BSMjMzPT6LwCAkpISMjMz6xxXXFxcY4w/x/d3v1AI9bkerV27dlx55ZVe/3ILt2Ce79tvv815551H165dA9ovnEJ9vkeL9PcbrHNNSEigf//+vPjii/Tq1Yu1a9eSlZUVtd9tbed7tKb+3a5evZp58+bx1FNPAeYf6YqKClq3bs3PP/8cdd9tfed7tEh/t7WJySsf2dnZ7Nq1i4qKCoDqy1jHXr46//zz2bJlS/XnrVu3Mnjw4HqP39D9QiHU53qsZs2accoppzR8wo0UrPPduHEjO3bs4NJLLw1ov3AL9fkeK5Lfb7D/WY6Li6Nv37507NgxoP3CJdTne6ym/N0OGDCATZs2sWHDBjZs2MDf/vY3+vfvz4YNG+rcL1JCfb7HivS/l32K9IrXSLniiiuMFStWGIZhGE8//bRxxx13GIZhGA888ED1CuENGzYYZ599tmEYhnH48GGjb9++RmFhoddx/vKXvxjjxo3z2ubPfuEUynN99dVXjU2bNhmGYRjffvutMW3atJCeiz8ae75bt241brnlFmP79u3G9u3bjQ0bNhiPPfZYvftFSijP12rfb2PPdePGjcb+/fsNwzCMAwcOGCNGjDA8Hk+9+0VKKM832r7boz399NNeT39E43d7tGPP12rfrS8xGz727t1r3HDDDcY999xj3HbbbUZ5eblhGIbRr18/Y8mSJdXjnnrqKeP22283br75ZuP999/3Osann35qnH322cbxxx9vrFq1yutnde0XbqE812nTphnJycnGxRdfbMyZM8dwu93hOak6NOZ8CwoKjE6dOhmA12vRokV17hdJoTxfq32/jf1n+brrrjPS0tKM66+/3pg2bZrx008/eR0/mr5bw6j7fKPtuz3asX+M/d0vnEJ5vlb7bn2xGYZhRO66i4iIiMSamFzzISIiIpGj8CEiIiJhpfAhIiIiYaXwISIiImGl8CEiIiJhpfAhIiIiYaXwISIiImGl8CEiIiJhpfAhIiIiYaXwISIiImGl8CEiIiJhpfAhIiIiYfX/ARZEPfPprcONAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "x = np.linspace(0,0.05,500)\n",
    "plt.plot(d,pc,'o')\n",
    "plt.plot(d,lognorm.cdf(d,params[0],scale=params[2]))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "14cf8ff8-26dc-488d-96cf-1e7d9cd1bc68",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "params = lognorm.fit(d,floc=0.0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "c38d847d-bc41-47a4-8fa6-fbddb7a21e1f",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(0.42293333848948156, 0.0, 0.021809608343306013)"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "params"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 62,
   "id": "59272be3-f545-42e5-b86e-e0709b1bacf9",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x149bfc79310>]"
      ]
     },
     "execution_count": 62,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAh8AAAGdCAYAAACyzRGfAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABLC0lEQVR4nO3deVxVdf7H8dcFgYsK1wUNKEfNMRuzQlLSxspMJ6zUaRMXzH3BJW0q06YyS1MnRzNLUdtcwWXMNc3MLUcFTakxS9OfpQKu2AVE1nt+fzgwoqBclnu58H4+Hvcx3ev3nvvhDHXffleTYRgGIiIiIg7i5uwCREREpHJR+BARERGHUvgQERERh1L4EBEREYdS+BARERGHUvgQERERh1L4EBEREYdS+BARERGHquLsAq5ls9lISEjAx8cHk8nk7HJERESkCAzDICUlhcDAQNzcbty3Ue7CR0JCAvXq1XN2GSIiIlIMJ0+e5Lbbbrthm3IXPnx8fIArxfv6+jq5GhERESmK5ORk6tWrl/c9fiPlLnzkDrX4+voqfIiIiLiYokyZ0IRTERERcSiFDxEREXEohQ8RERFxKIUPERERcSiFDxEREXEohQ8RERFxKIUPERERcSiFDxEREXEohQ8RERFxqGKFjw0bNhASEsKvv/5aaJs1a9YwcuRIhgwZwubNm4tbn4iIiFQwdm+vfvbsWbKzs9m7d2+hbX7++WcmTJhATEwMhmHQokUL1q5dy6233lqiYkVERMT12d3zUbduXZ544okbtnn//fcJDQ3FZDLh5uZG69atmT17drGLFBERkYqjWMMubm43ftuWLVuoX79+3vPGjRuzffv2AttmZGSQnJyc7yEiIiKlLykpiaeffppvvvnGqXWUyYTT+Ph4atWqlffcx8eHhISEAttOmjQJi8WS96hXr15ZlCQiIlKp7d69m+bNm/PFF1/Qv39/srKynFZLmYQPk8mE2WzOe56ZmYmHh0eBbceOHYvVas17nDx5sixKEhERqZRsNhvvvfceDz30ECdOnOCPf/wjX3zxRaHfy45g94TToggMDMRqteY9T0lJITAwsMC2Xl5eeHl5lUUZIiIildr58+fp3bs3X375JQDdunVjzpw5+Pr6OrWuMgkfjz76KEeOHMl7fvToUR555JGy+CgREREpwLfffkv37t2Jj4/HbDbzwQcfMGDAAEwmk7NLK96wi2EY+f4XYOrUqRw+fBgg394e2dnZxMbGMnDgwJLWKiIiIjdhs9mYOHEibdu2JT4+niZNmhATE8PAgQPLRfCAYvR8pKamsnDhQgDmz5/P8OHD8fPzIyoqigYNGtCkSRPuvfde+vbty8svv0xmZibTp0/H39+/1IsXERGR/zlz5gy9evXi66+/BqBXr17MmjWL6tWrO7my/EzG1d0X5UBycjIWiwWr1er0MSkRERFXsXXrVnr06MHp06fx9vbmo48+ok+fPg7r7bDn+1tnu4iIiLiwnJwcxo8fT/v27Tl9+jRNmzZl37599O3bt9wMs1yrTCacioiISNlLTEykZ8+ebN26FYB+/foxc+ZMqlat6uTKbkzhQ0RExAV9/fXXhIeHc/bsWapVq0ZkZCTh4eE3fE+OzSD2eBJnU9Kp62MmpGEt3N0c3zui8CEiIuJCsrOzeeutt3j33XcxDIN77rmHpUuXcuedd97wfRsPJjJ+7SESrel5rwVYzIzr1JTQZgFlXXY+mvMhIiLiIk6dOkW7du2YOHEihmEwePBg9uzZU6TgEbFof77gAXDamk7Eov1sPJhYlmVfR+FDRETEBWzYsIGgoCC+/fZbfHx8iI6OJjIyEm9v7xu+L8dmMH7tIQpa2pr72vi1h8ixOW7xq8KHiIhIOZaVlcWrr77K448/zoULFwgODmb//v2EhYUV6f2xx5Ou6/G4mgEkWtOJPZ5UShXfnOZ8iIiIlFMnTpygW7du7N69G4Dhw4czdepUu85EO5tSePAoTrvSoPAhIiJSDq1Zs4Y+ffpw8eJFLBYLn3zyCc8884zd16nrY755IzvalQYNu4iIiJQjmZmZ/O1vf6NLly5cvHiRli1bcuDAgWIFD4CQhrUIsJgpbEGtiSurXkIa1ip2zfZS+BARESknjh8/Tps2bZg+fToAL774Ijt37qRhw4bFvqa7m4lxnZoCXBdAcp+P69TUoft9KHyIiIiUAytXrqR58+bs3buXmjVrsnr1aqZNm4anp2eJrx3aLIDZ4cH4W/IPrfhbzMwOD3b4Ph+a8yEiIuJEGRkZvPzyy3z44YcAtG7dmqioKOrXr1+qnxPaLIAOTf21w6mIiEhldvToUcLCwti/fz8Ao0ePZsKECXh4eBTYvqTbo7u7mWjdqHap1F4SCh8iIiJOsHTpUgYOHEhKSgq1a9dmwYIFPP7444W2L0/bo5eU5nyIiIg40OXLlxkyZAjdunUjJSWFBx98kLi4uJsGj/K0PXpJKXyIiIg4yOHDh2nVqhVz5szBZDLRe+jf+HrzN9x2222Fvqc8bo9eUgofIiIiDrBo0SKCmgfzww8/4Fa1BnWee5ttPu1o+88dN+y5KI/bo5eUwoeIiEgZSktLo3///vTq1Yv0y2l4/eEeAvp+gHfD5sDNh07K4/boJaUJpyIiImXk0KFDPPfccxw6dAhMJiwPdMfyQBgmN/e8NgZXNvsav/YQHZr6X7d6pTxuj15S6vkQEREpZYZh8Nlnn9GiRQsOHTpE7Tp1uSVsIjXa9MgXPPLaU/jQSXncHr2kFD5ERERKUWpqKr1796Zfv35cvnyZDh068M+orzDXv+em7y1o6KQ8bo9eUgofIiIidsixGew+doHVcfHsPnYh3yqTH374gZYtW7Jw4ULc3NyYOHEiGzdu5I4G9Yp07cKGTsrb9uglpTkfIiIiRVTYRl9vPvknTu5ex8iRI0lPT+fWW28lKiqKBx98EPjf0Mlpa3qBS2ZNXAkSNxo6KU/bo5eUwoeIiEgR5G70dW14SDibxHNh3bn00w4AOnbsyIIFC/Dz88trkzt0ErFoPybIdw17hk7Ky/boJaVhFxERkZsobKOvzDPHSJg/8krwcHNj8pQprFu3Ll/wyFXRhk5KQj0fIiIiN3HtRl+GYZB6YD1JWz6GnGzcfetQp/NoHnqmL25uhf+9viINnZSEwoeIiMhNXL0KxZaeyoWNM0k7/G8AvP94P7UfH4W7t0+RNvqqKEMnJaHwISIilUpxjqXPXYWSkXiE86unkG09A25VqNm2Lz4tOmMymfK1kxtT+BARkUqjuMfSt2xQEw6u5/SGeWDLxt1yC3W6vIpXwB1A0VaryP9owqmIiFQKxT2WPikpiWefeZrf1s8GWzZV73iAwD4z8gUPcL2NvpxJPR8iIuJy7B06udmx9IWdrbJnzx7CwsI4ceIEnp6eDHh5HN9VC+F0ckZeG/8i9JxIfgofIiLiUoozdGLPsfStG9XGZrPxz3/+k9dee43s7GwaNWrEsmXLCA4OLtacEclP4UNERFxGYRt95Q6dFLZfhj3H0p8/f54+ffqwfv16AMLCwpg7dy6+vr6AVquUBs35EBERl3CzoRO4MnRy9VkruYq6CiXxcBxBQUGsX78eLy8v5syZQ1RUVF7wkNKh8CEiIi7BnqGTa93sWHoMG7YDKxneswvx8fE0adKE2NhYBg0alLeMVkqPwoeIiLgEe4ZOrnWjY+ltl37nzLJxnNz0KTk5OYSHh7Nv3z7uueeekpYshVD4EBERl1DUoRN7jqVPP/EDp+ePJP3XA3h7e/Ppp5+yYMECqlevXio1S8E04VRERFxCaR5Lv/voOT6aPoVlS9/HZrPRtGlTli1bxl133VVm9cv/qOdDRERcwo2GTuzZ6Ovc2TOMi+hOdOQ0bDYbffv2JTY2VsHDgRQ+RETEZZT0WPrNmzdz7733smXLFqpVq8aCBQv49NNPqVatWlmWLdfQsIuIiDhcSTbqKs6x9NnZ2bz11lu8++67GIbB3XffzbJly7jzzjtL60cSOyh8iIiIQxX3cLer2bPRV3x8PD169GDHjh0ADB48mOnTp+Pt7W1/8VIqNOwiIiIOU9zD3Yr9eRs3EhQUxI4dO/Dx8SEqKorIyEgFDydT+BAREYcoyQ6l9srKymLMmDF07NiR8+fP07x5c7777ju6detW4mtLySl8iIiIQ5Rkh1J7nDhxgrZt2zJlyhQAhg0bxq5du2jcuHGJriulR3M+RETEIUqyQ2lRrV27lt69e3Px4kUsFguffPIJzzzzTLGvJ2VDPR8iIuIQJd2h9EYyMzN56aWX6Ny5MxcvXqRly5bs379fwaOcUvgQERGHuNnhbiaurHq50Q6lBTl+/DgPPvgg06ZNA2DUqFHs3LmT22+/vWQFS5lR+BAREYcorR1Kr7Zy5UqaN29ObGwsNWrUYNWqVUyfPh1PT8/SKVrKhMKHiIg4TEl3KM2VkZHBiBEjeOaZZ7BarbRq1Yq4uDi6dOlSFmVLKdOEUxERcaji7FB6taNHjxIWFsb+/fsBeOWVV5g4cSIeHh5lWbaUIoUPERFxOHt2KL3asmXLGDBgACkpKdSuXZv58+fzxBNPlEGFUpY07CIiIuXe5cuXiYiIICwsjJSUFNq0aUNcXJyCh4tS+BARkXLt8OHDtGrVisjISEwmE6+99hpbt27ltttuc3ZpUkwadhERkXJr8eLFDB48mEuXLlGnTh0WLVrEX/7yF2eXJSWkng8RESl30tLSGDBgAOHh4Vy6dIm2bdvy/fffK3hUEHb3fFy6dInRo0dTs2ZNUlNTmTJlCl5eXte1e/fdd7HZbLi5uZGWlsY777yDyVT0tdsiIlK+5diMYq9YuZFDhw7RtWtXfvzxR0wmE2+++SZvvPEG7u7upVC1lAd2h4+IiAieeuopnnrqKRYsWMDYsWPzdpXLtW7dOo4fP868efMA6NevH8uWLSMsLKx0qhYREafaeDCR8WsP5TsoLsBiZlynpkXeq6Mgn3/+OcOGDSMtLQ1/f38WL15Mu3btSqNkKUfsGnZJSEhg+fLldOzYEYCOHTsSGRlJSkpKvnY//fRTvte8vb2xWq2lUK6IiDjbxoOJRCzaf90Jtaet6UQs2s/Gg4l2XzM1NZXevXvTt29f0tLSaN++PXFxcQoeFZRd4WPbtm34+flhNl/Zma5OnTp4enoSGxubr13nzp1Zs2YNK1as4OLFi5w/f55evXoVeM2MjAySk5PzPUREpHzKsRmMX3sIo4A/y31t/NpD5NgKalGw//znP7Rs2ZIFCxbg5ubGhAkT+Oqrr7jllltKpWYpf+wKH/Hx8dSqlf/AHx8fHxISEvK91qRJE5YsWcLzzz9Pjx49+Oyzz/D29i7wmpMmTcJiseQ96tWrZ+ePICIixZFjM9h97AKr4+LZfexCkQJD7PGk63o8rmYAidZ0Yo8n3fRahmEwb948QkJC+PnnnwkMDGTr1q38/e9/x81N6yEqMrv+3zWZTHm9HrkyMzML3NI2PT2dNWvWcOzYMbp37052dnaB1xw7dixWqzXvcfLkSXtKEhGRYth4MJE2U7bQfd4eRkbH0X3eHtpM2XLTIZOzKYUHD3vapaSk0LNnTwYNGkR6ejodO3YkLi6Ohx56qMg/g7guu8JHYGDgdXM3UlNTCQwMzPfatm3b2LdvH+3bt2fbtm0cOHDgukmpuby8vPD19c33EBGRslOSORt1fcyF/llR2x04cIDg4GCioqJwd3dnypQprFu3jjp16hTtBxCXZ1f4aNu2LadOnSIzMxMgb7glJCQkX7sVK1bQuHFj4EpgmTBhAjt27CiNekVEpARKOmcjpGEtAixmCltQa+LKqpeQhrWu+zPDMJg1axatW7fm6NGj1KtXjx07djB69GgNs1Qydvd8hIaGsn37dgA2bdrE0KFDMZvNTJ06lcOHDwMQFBTEgQMH8t5nMpmuCygiIuJ4JZ2z4e5mYlynpgDXBZDc5+M6Nb1uvw+r1UrXrl0ZNmwYGRkZdOrUibi4OB544IFi/iTiyuze5yMyMpIxY8YQExNDUlISkydPBiAqKooGDRrQpEkT+vXrx4kTJ5g8eTJ+fn6cPHmS1157rdSLFxER+5TGnI3QZgHMDg++bp8P/0L2+di3bx9du3bl+PHjeHh4MGXKFEaNGqWNJysxk2EYRV8P5QDJyclYLBasVqvmf4iIlLLdxy7Qfd6em7aLGtjqpkfe32yHU8Mw+OCDD3jllVfIysqiQYMGLF26VD3hFZQ93986WE5EpBLJnbNx2ppe4LwPE1d6MAqas3EtdzdToQHl4sWL9OvXj1WrVgHw9NNP88knn1CjRo1i1y4Vh2b4iIhUIsWds2GPPXv20Lx5c1atWoWnpyczZ85kxYoVCh6SR+FDRMSFFWejsNw5G/6W/Mth/S1mZocHF/tsFpvNxtSpU3nwwQf57bffaNSoEbt27WL48OGa3yH5aNhFRMRFleRwt9BmAXRo6l9qp9JeuHCB3r17s379egDCwsKYO3eu5u5JgTThVETEBeVuFHbtf8Bzo0NJejDstXPnTrp3786pU6fw8vJixowZDBo0SL0dlYw9398adhERcTFlcbhbcdhsNiZNmpS3AeUdd9xBTEwMgwcPVvCQG1L4EBFxMaV5uFtxnT17lscff5zXXnuNnJwcwsPD+e6777j33nvL7DOl4tCcDxERF1Nah7sV1/bt2+nevTuJiYl4e3vz4Ycf0rdvX/V2SJGp50NExMWUxuFuxZGTk8Pbb79Nu3btSExM5E9/+hN79+6lX79+Ch5iF/V8iIi4mNLcKKyoTp8+Tc+ePdmyZQsAffv2ZebMmVSrVq3UPkMqD/V8iIg4UXH26XDERmFX27x5M0FBQWzZsoWqVauyYMECPv30UwUPKTb1fIiIOElJ9+mw53C34sjOzmb8+PFMnDgRwzC4++67WbZsGXfeeWeJry2Vm/b5EBFxgtLap+Nmh7sVV3x8PD169GDHjh0ADBw4kBkzZuDt7V3ia0vFpIPlRETKsZvt02Hiyj4dHZr63zRI3Ohwt+LauHEjvXr14vz581SvXp25c+fSvXv3Uv0Mqdw050NExMHKwz4dBcnKymLs2LF07NiR8+fPExQUxP79+xU8pNSp50NExMGcvU9HQU6ePEm3bt3YtWsXAMOGDWPq1KmYzaW7XFcEFD5ERBzOWft0FGbdunX07t2bpKQkfH19+eSTT3j22Wcd8tlSOWnYRUTEwXL36ShsNoeJK6teSnOfjoJkZmby0ksv0alTJ5KSkmjRogUHDhxQ8JAyp/AhIuJgjt6noyC//vorDz74INOmTQNg1KhR7Ny5k9tvv73MPlMkl8KHiIgT5O7T4W/JP7TibzEXeZltcX3xxRc0b96c2NhYatSowapVq5g+fTpeXl5l9pkiV9OcDxERJwltFkCHpv5lsk9HQTIyMnjllVeYOXMmAK1atSI6Opr69euXyeeJFEbhQ0TEicpin46CHDt2jLCwML777jsAXn75Zd599108PDzK/LNFrqXwISJSwS1fvpwBAwaQnJxM7dq1mT9/Pk888YSzy5JKTHM+REQqqPT0dIYOHUrXrl1JTk6mTZs2xMXFKXiI0yl8iIhUQEeOHKFVq1bMnj0bgLFjx7J161Zuu+02J1cmomEXEZEKZ8mSJQwePJjU1FTq1KnDwoULeeyxx5xdlkge9XyIiFQQaWlpDBw4kJ49e5Kamkrbtm2Ji4tT8JByR+FDRKQC+Omnn7j//vv5+OOPMZlMvPnmm2zevJnAwEBnlyZyHQ27iIi4uPnz5zN06FDS0tK45ZZbWLJkCe3atXN2WSKFUs+HiIiLunTpEr1796ZPnz6kpaXRvn17vv/+ewUPKfcUPkREXNB//vMfWrRowYIFC3Bzc+Odd95h48aN3HLLLc4uTeSmNOwiIuJCDMPgk08+YcSIEaSnpxMYGMiSJUt4+OGHnV2aSJEpfIiIuIiUlBSGDBnCkiVLAAgNDWXBggXUqVPHyZWJ2EfDLiIiLiAuLo777ruPJUuW4O7uzuTJk1m/fr2Ch7gk9XyIiJRjhmEQGRnJiy++SEZGBvXq1SM6OpoHHnjA2aWJFJvCh4hIOWW1Whk4cCDLly8HoFOnTnz22WfUrl32p+CKlCUNu4iIlEP79u0jODiY5cuXU6VKFf75z3+yevVqBQ+pENTzISJSjhiGwcyZM3n55ZfJysqifv36LF26lPvvv9/ZpYmUGoUPEZFy4uLFi/Tr149Vq1YB8NRTT/HJJ59Qs2ZN5xYmUso07CIiUg7ExMTQvHlzVq1ahaenJx988AH/+te/FDykQlL4EBFxIsMw+Oc//0mbNm347bffuP3229m1axcjRozAZDI5uzyRMqFhFxERJ7lw4QJ9+vRh3bp1AHTt2pW5c+disVicXJlI2VLPh4iIE/z73/8mKCiIdevW4eXlxezZs4mOjlbwkEpB4UNExIFsNhuTJ0/m4Ycf5tSpU9xxxx3ExMQwZMgQDbNIpaFhFxERBzl37hzPP/88GzduBKBnz57Mnj0bHx8fJ1cm4lgKHyIiDrB9+3Z69OhBQkIC3t7ezJw5k379+qm3QyolDbuIiJShnJwc3nnnHdq1a0dCQgJ/+tOfiI2NpX///goeUmmp50NEpIycPn2a8PBwvvnmGwD69OnDhx9+SLVq1ZxcmYhzKXyIiJSBb775hp49e3LmzBmqVq3K7Nmzef75551dlki5oGEXEZFSlJOTw5tvvkmHDh04c+YMzZo1Y9++fQoeIldRz4eISClJSEigR48ebN++HYCBAwcyY8YMvL29nVyZSPmi8CEiUgq++uorwsPDOX/+PNWrV2fOnDn06NHD2WWJlEsadhERKYHs7GxeHTOG0NBQzp8/T+M/3U3s3n0KHiI3oPAhIlJMJ0+e5N6QB/jHlCkAVG/+BBmPj6ffF6fYeDDRydWJlF8KHyIixbB+/XruuvteDh3Yi8mzKn5dxlD7LxGYqnhy2ppOxKL9CiAihVD4EBGxQ1ZWFi+//DJPPvkkKdaLePr/kYA+M6h2Z5u8NsZ//3f82kPk2IyCLyRSidk94fTSpUuMHj2amjVrkpqaypQpU/Dy8iqwbVpaGnPmzKF27drccccdtGrVqsQFi4g4y6+//kq3bt2IiYkBwOe+ztRs2xdTFY/r2hpAojWd2ONJtG5U28GVipRvdvd8RERE0L59eyZMmEBwcDBjx44tsF1SUhLPPfccf/3rX3n++ecVPETEpa1atYrmzZsTExNDjRo1GDPtY2q1H1Rg8Lja2ZR0B1Uo4jrsCh8JCQksX76cjh07AtCxY0ciIyNJSUm5rm23bt145ZVXaNiwYelUKiLiBBkZGYwcOZKnnnqK33//nfvvv58DBw7QufNfi/T+uj7msi1QxAXZFT62bduGn58fZvOVf5nq1KmDp6cnsbGx+dqtW7eOX375hdjYWB5//HHGjh1LVlZWgdfMyMggOTk530NEpDw4duwYf/7zn/nggw8AeOmll9ixYwcNGjQgpGEtAixmCjsazgQEWMyENKzlsHpFXIVd4SM+Pp5atfL/i+Tj40NCQkK+1xYvXsz999/PqFGjWLhwIYsXL+aNN94o8JqTJk3CYrHkPerVq2fnjyAiUvqWL19OcHAw3333HbVq1WLt2rVMnToVT09PANzdTIzr1BTgugCS+3xcp6a4u+nkWpFr2RU+TCZTXq9HrszMTDw88o95/vjjj7Rp0wZPT09q167NwIEDWbBgQYHXHDt2LFarNe9x8uRJO38EEZHSk56eztChQ+natSvJycn8+c9/Ji4ujieffPK6tqHNApgdHoy/Jf9/F/0tZmaHBxPaLMBRZYu4FLtWuwQGBmK1WvO9lpqaSmBgYL7XsrOzycnJyXt+zz33kJSUVOA1vby8Cl0tIyLiSL/88gtdu3YlLi4OuPKXo/Hjx1/3F6yrhTYLoENTf2KPJ3E2JZ26PleGWtTjIVI4u8JH27ZtGTRoEJmZmXh6euYNt4SEhORrd8899/DLL7/870OqVOHOO+8shXJFRMpGVFQUgwYNIjU1lTp16rBw4UIee+yxIr3X3c2k5bQidrBr2CUwMJDQ0NC8Exs3bdrE0KFDMZvNTJ06lcOHDwMwcuRI1q9fT0ZGBgC7du1ixIgRpVy6iEjJXb58mYEDB9KjRw9SU1N5+OGHiYuLK3LwEBH72b3JWGRkJGPGjCEmJoakpCQmT54MXPlbQ4MGDWjSpAmtW7fmnXfeYcSIETRp0gSz2Uz//v1LvXgRkZL46aef6Nq1KwcPHsRkMvHGG2/wxhtvUKWKDvwWKUsmwzDK1d6/ycnJWCwWrFYrvr6+zi5HRCqoBQsWEBERQVpaGrfccguLFy/m0UcfdXZZIi7Lnu9vne0iIpXKpUuX6Nu3L7179yYtLY1HH32UuLg4BQ8RB1L4EJFK4+DBg7Rs2ZLPP/8cNzc33n77bb766iv8/f2dXZpIpaKBTRGp8AzD4NNPP2X48OGkp6cTGBjIkiVLePjhh51dmkilpPAhIhVaSkoKERERLF68GIDHHnuMhQsXUqdOHSdXJlJ5adhFRCqs77//nhYtWrB48WLc3d2ZNGkSX375pYKHiJOp50NEKhzDMJgzZw6jRo0iIyOD2267jejoaP785z87uzQRQeFDRCoYq9XKoEGDWLZsGQBPPvkkn3/+ObVrF74DaY7N0PboIg6k8CEiFcZ3331HWFgYx44do0qVKkyZMoUXX3wRk6nwILHxYCLj1x4i0Zqe91qAxcy4Tk11MJxIGdGcDxFxeYZhMHPmTB544AGOHTtG/fr12blzJ3/7299uGjwiFu3PFzwATlvTiVi0n40HE8u6dJFKSeFDRFzaxYsXeeaZZ3jhhRfIzMzkr3/9KwcOHOD++++/4ftybAbj1x6ioC2ec18bv/YQObZytQm0SIWg8CEiLismJobg4GC++OILPDw8mDFjBitXrqRmzZo3fW/s8aTrejyuZgCJ1nRijyeVYsUiAgofIuKCDMNg2rRptGnThl9//ZXbb7+dXbt28cILL9xwmOVqZ1MKDx7FaSciRacJpyLiUi5cuECfPn1Yt24dAM899xzz5s3DYrHYdZ26PuZSbSciRaeeDxFxGbt27aJ58+asW7cOLy8vZs2axdKlS+0OHgAhDWsRYDFTWD+JiSurXkIa1ipRzSJyPYUPESn3bDYbU6ZM4aGHHuLkyZM0btyYPXv2EBERUeRhlmu5u5kY16kpwHUBJPf5uE5Ntd+HSBlQ+BCRcu3cuXM88cQTjBkzhpycHHr06MF3331HUFBQia8d2iyA2eHB+FvyD634W8zMDg/WPh8iZURzPkSk3NqxYwfdu3cnISEBs9nMhx9+SL9+/Yrd21GQ0GYBdGjqrx1ORRxI4UNEyp2cnBwmTZrEuHHjsNls3HnnnSxfvpxmzZqVyee5u5lo3ajw7ddFpHQpfIhIuXLmzBnCw8PZvHkzAL179+ajjz6iWrVqTq5MREqLwoeIlBtbtmyhR48enDlzhqpVqzJr1ix69+5dpPfqcDgR16HwISJOl5OTw9tvv80777yDYRg0a9aMpUuX0rRp0yK9X4fDibgWrXYREadKSEigffv2vP322xiGwYABA4iJibEreOhwOBHXovAhIk7z1VdfERQUxLZt26hevTqLFy9m3rx5VK1atUjv1+FwIq5J4UNEHC47O5vXXnuN0NBQzp07x7333st3331Hjx497LqODocTcU2a8yEiDnXq1Cm6d+/Ozp07AYiIiGDatGmYzfafoaLD4URck8KHiDjM+vXr6d27NxcuXMDX15d58+bRtWvXYl9Ph8OJuCYNu4hImcvKyuKVV17hySef5MKFC9x3333s37+/RMEDdDiciKtS+BCRMvXbb7/x0EMPMXXqVABeeOEF/v3vf9OoUaMSX1uHw4m4JoUPESkzq1evJigoiD179lCjRg1WrlzJjBkz8PLyKrXP0OFwIq5Hcz5EpNRlZmYyevRoZsyYAUBISAhLly6lQYMGZfJ5OhxOxLUofIhIqfq///s/wsLC2LdvHwAvvfQS7777Lp6enmX6uTocTsR1KHyISKlZsWIF/fv3Jzk5mVq1avH555/TqVOnIr1XZ7OIVB4KHyJSYunp6bz00kvMmjULgAceeIDo6Gjq1atXpPfrbBaRykUTTkWkRH755Rdat26dFzzGjBnDtm3b7AoeOptFpHJR+BCRYouKiiI4OJi4uDj8/PzYsGEDkyZNwsPDo0jv19ksIpWTwoeI2O3y5csMGjSIHj16kJqaykMPPURcXByhoaF2XUdns4hUTprzISJ2+fnnn+natSv/+c9/MJlMvP7667z55puY3NzZfeyCXRNGdTaLSOWk8CEiRbZgwQIiIiJIS0vjlltuYdGiRbRv377YE0Z1NotI5aRhFxG5qUuXLtG3b1969+5NWloa7dq1Iy4uLi94FHfCqM5mEamcFD5EJE+OzWD3sQusjotn97EL5NgMfvzxR0JCQvj8889xc3Pj7bffZtOmTfj7+5d4wqjOZhGpnDTsIiLA9XttGIZBlaPbiP/yIzLS0wkICGDJkiW0bds27z32TBgtbPfR3LNZrh228dc+HyIVlsKHSAVS3F1Cc4dOcvsnbBlpJG2axaVD2wAIfuBhNnyxjLp16+Z7X2lNGNXZLCKVi8KHSAVR3Emf1w6dZJ79P86tnkJ2UjyY3KjxUC88OvSktl+d695bmhNGdTaLSOWhOR8iFUBJJn3mDp0YhkFK3AYSF7xEdlI87j5+3NJjEpZWz3E6JbPAvTY0YVREikPhQ8TFlXTS59mUdGwZaZxf8w+SvvoIcrLwbtSSgL4fYL7trnztrqUJoyJSHAofIi6upLuEXvztMImfjyTt52/BzZ2aj/SjzjNv4O7tm69dYUMnuRNG/S35/9zfYmZ2eLAmjIrIdTTnQ8TFFXfSp2EYfPTRR7z00ktkZ2bi7luXOp1H43XrnfnambgSJG40dKIJoyJiD4UPERdXnEmfv//+O/3792flypUAtG4Xyqm7++Burp5v+MaeoRNNGBWRotKwi4iLs3fSZ2xsLM2bN2flypV4eHgwY8YM/r35S+YOeEhDJyLiEOr5EHFxuZM+IxbtxwSF9ly4mWD69Om8+uqrZGVlcfvtt7N06VJatGgBaOhERBxH4UOkArjZLqEhgV506dKFtWvXAvDss8/y8ccfY7FY8l1HQyci4ggKHyIVRGE9FzF7dhP0eDdOnjyJl5cX06dPZ8iQIZhM6tEQEedQ+BCpQK7uubDZbEyd+h6vvfYaOTk5NG7cmGXLlhEUFOTcIkWk0lP4EKmAzp07R+/evdmwYQMA3bt3Z86cOfj4+Di5MhERhQ+RCmfHjh10796dhIQEzGYzM2fOpH///hpmEZFyQ0ttRSoIm83GxIkTeeSRR0hISODOO+8kNjaWAQMGKHiISLlid/i4dOkSw4YN4/XXX2fUqFFkZGTcsH10dDRt27Ytbn0iUgRnzpwhNDSU119/HZvNxvPPP8/evXu5++67nV2aiMh17A4fERERtG/fngkTJhAcHMzYsWMLbZuQkMDEiRNLVKCI3NiWLVsICgri66+/pmrVqnz22WfMnz+f6tWrO7s0EZEC2RU+EhISWL58OR07dgSgY8eOREZGkpKSUmD7yZMnExERUfIqReQ6OTk5vPXWW7Rv357Tp09z1113sXfvXvr06ePs0kREbsiu8LFt2zb8/Pwwm69swVynTh08PT2JjY29ru28efMIDw+natWqN7xmRkYGycnJ+R4icmMJCQm0b9+e8ePHYxgG/fv3JzY2lqZNmzq7NBGRm7IrfMTHx1OrVv6TLX18fEhISMj32tGjR0lOTiYkJOSm15w0aRIWiyXvUa9ePXtKEql0Nm3aRFBQENu2baNatWosWrSIjz/++KZBX0SkvLArfJhMprxej1yZmZl4eHjkPc/JyWHWrFmMGjWqSNccO3YsVqs173Hy5El7ShKpcHJsBruPXWB1XDy7j10gx3bltJbs7Gz+/ve/Exoayrlz57j33nvZv38/PXv2dHLFIiL2sWufj8DAQKxWa77XUlNTCQwMzHu+a9cuIiMj+fTTT4Er4SQzM5MaNWrw+++/X3dNLy8vvLy8ilG6SMWz8WDideezBFjMDG1Zg7njR7Fz504AhgwZwrRp0/D29nZWqSIixWYyDMO4ebMrEhISuOOOO0hKSsLT05OEhAQaN27MhQsX8npE0tPTOX36dN57VqxYwYoVK4iOjqZBgwY3/Yzk5GQsFgtWqxVfX1/7fyIRF7XxYCIRi/Zz7b+Ql4/t5fz66dguJ+Pj48PHH39M165dnVKjiEhh7Pn+trvnIzQ0lO3bt9OhQwc2bdrE0KFDMZvNTJ06lU6dOtGkSZN8ISN3gmpRgodIZZVjMxi/9lC+4GHkZPP7jgUkx64EoFpgY/Zt+5I7Gv/ROUWKiJQSu7dXj4yMZMyYMcTExJCUlMTkyZMBiIqKokGDBjRp0qTUixSp6GKPJ+Ubasm2nuXcmilkJhwGwOe+TtRs248LbjWdVaKISKmxa9jFETTsIpXR6rh4RkbHAZD2yx4ufPk+tvRUTF7V8Os4kqpNHgBgRrcgugTd6sRKRUQKVmbDLiJSNur6mDFysri47XNS9q0GwDPgDvw6j8ajhn++diIirk7hQ6QcqIOVC9FjuHTqv8MsLf9KzYd7Y3K/sozdBPhbzIQ0rHWDq4iIuAaFDxEn+9e//kX//v25ZLXiZq6O3xMv4v3H+/P+PPc82nGdmuLuptNpRcT12X2wnIiUjvT0dIYPH86zzz6L1WrlgQce4LM1W7j9vofztfO3mJkdHkxoswAnVSoiUrrU8yHiBL/88gthYWEcOHAAgFdffZV33nkHDw8Pej5qEHs8ibMp6dT1uTLUoh4PEalIFD5EHCw6OppBgwaRkpKCn58fCxcuJDQ0NO/P3d1MtG5U24kVioiULQ27iDjI5cuXGTx4MN27dyclJYWHHnqIuLi4fMFDRKQyUPgQcYCff/6Z+++/n7lz52IymXj99df55ptvuPVW7dkhIpWPhl1EytjChQuJiIjg0qVL1K1bl8WLF9O+fXtnlyUi4jTq+RApI5cuXaJfv348//zzXLp0iXbt2hEXF6fgISKVnsKHSBn48ccfCQkJ4bPPPsPNzY3x48ezadMmAgK0XFZERMMuIqXIMAw+/fQzhg0fTkb6Zfzq3kJ0VBSPtnvE2aWJiJQbCh8ipSQ1NZUu3fuwZd2/ADA3aI75yZd4PcZGVt1EbRImIvJfGnYRKUCOzWD3sQusjotn97EL5NhufPjzDz/8wJ1333sleJjcqPHQ89TtOh73ajU4bU0nYtF+Nh5MdFD1IiLlm3o+RK6x8WAi49ceItGanvdagMXMuE5Nr+u9MAyDuXPnMnLkSDIyMnCvXhu/LqMx33bX/9pw5XyW8WsP0aGpv3YrFZFKTz0fIlfZeDCRiEX78wUPoMDei+TkZLp3786QIUPIyMjA+/YWBPT9IF/wyGUAidZ0Yo8nlfWPICJS7qnnQ+S/cmwG49ceoqABlmt7L76PO0BYWBhHjx6lSpUqhI8YwxaPEEymG+f5synpN/xzEZHKQD0fIv8Vezzpuh6PqxlAwu+XGT3+H7Ru3ZqjR4/yhz/8gW+//ZZBw0bdNHgA1PUxl2LFIiKuST0fIv91s14JW3oqFzZ8wLQjuwDo0qULn376KbVq1SLHZhBgMXPaml5gz4kJ8LdcOaFWRKSyU8+HyH/dqFciI/EIiZ+PJO3ILqp4ePD+++/zxRdfUKvWlTDh7mZiXKemwJWgcbXc5+M6NdVkUxERFD5E8oQ0rEWAxZwvPBiGQfLe1ZxeNJps6xm8avrz7bc7GTlyJCZT/iAR2iyA2eHB+Fvyhxh/i5nZ4cHa50NE5L807CLyX7m9FxGL9mMCsi+ncOHL97l8NAaAqnc8wPzPPqHV/XcWeo3QZgF0aOpP7PEkzqakU9fnylCLejxERP7HZBjGjXdPcrDk5GQsFgtWqxVfX19nlyOV0MaDibz80Qp+XvIOOcnnwL0KDZ6IYNaEMXS8O9DZ5YmIlEv2fH+r50PkKjabjR++XMjP814jJyeHgHoNePfDT+n1ZFv1XoiIlBKFD5H/On/+PM8//zwbNmwAoFu3bsyZM0c9cCIipUwTTkWAb7/9lqCgIDZs2IDZbGbu3LksWbJEwUNEpAwofEilZrPZmDhxIm3btiU+Pp4mTZoQExPDwIEDr1vNIiIipUPDLlJpnTlzhl69evH1118D0KtXL2bNmkX16tWdXJmISMWm8CGV0tatW+nRowenT5/G29ubWbNm0adPH2eXJSJSKWjYRSqVnJwcxo8fT/v27Tl9+jR33XUX+/btU/AQEXEg9XxIpZGYmEjPnj3ZunUrAP369WPmzJlUrVrVyZWJiFQuCh9SKXz99deEh4dz9uxZqlWrRmRkJOHh4c4uS0SkUlL4kAopx2YQezyJxIuprPnsAxbMno5hGNxzzz0sW7aMJk2aOLtEEZFKS+FDKpyNBxMZv/YQJ0+e4vza98g49SMAjz/XixXz5+Dt7e3kCkVEKjeFD6lQNh5MJGLRftKO7eX8+unYLidj8vTGL3QEh25/iO3Hfie0mcKHiIgzKXxIhZFjMxi36geStn5KcuxKADxvaYRfl1fxqHnlQLjxaw/Roam/zmkREXEihQ+pMNb8+we+nz2KjISfAfC5rxM12/bDVMUDAANItKYTezyJ1o1qO7FSEZHKTeFDKoQ1a9bwfK/nyUi2YvKqhl/HkVRt8kCBbc+mpDu4OhERuZo2GROXlpmZyd/+9je6dOlCarIVz4DGBPSZUWjwAKjrY3ZghSIici31fEi5lbtc9mxKOnV9zIQ0rJVvrsbx48cJCwtj7969AIwcNYrdtR7j7KUcjAKuZwL8LVeuIyIizqPwIeVS7nLZROv/hkgCLGbGdWpKaLMAVq5cSb9+/bBardSsWZPPP/+czp075612MUG+AJIbWcZ1aqrJpiIiTmYyDKOgvyQ6TXJyMhaLBavViq+vr7PLESfIDRDX/mKaACM7i3sS1rA26jMAWrduTXR0NH/4wx/yvf9GwUVEREqfPd/f6vmQciXHZjB+7aECh00yLyZwfvUUfjtzDIDRo0czYcIEPDw88rULbRZAh6b+NxyyERER51H4kHIl9nhSvh6LXJd+2sGFjTMxMi/j5u3LP2bO4aX+3Qq9jrubSctpRUTKKYUPKVeuXQZry8rg4pZ5pMZtBMDrtrvw6/QKf7zvQWeUJyIipUDhQ8qVq5fBZl04xbnVk8k69ytgwtK6K5Y2PTC5uWu5rIiIC1P4kHIlpGEtAixmju7awIWvPsLISsetag38nnwJ74bNtVxWRKQCUPiQciUj/TK+ez/m/LpoALz+cA9+nV6mSvVaWi4rIlJBKHxIuXHo0CGee+45Dh06hMlk4tZ2vXALfgaTmztwpcdDy2VFRFyfwoc4nWEYfP755wwbNozLly/j7+/PkiVLeOjhtlouKyJSASl8iFOlpqYydOhQFi5cCECHDh1YtGgRdevWBdByWRGRCkgHy4nT/PDDD7Rs2ZKFCxfi5ubGxIkT2bhxY17wEBGRikk9H+JwhmEwb948Ro4cSXp6OrfeeitRUVE8+KD27hARqQwUPsShkpOTGTx4MNHRV1azdOzYkQULFuDn5+fkykRExFE07CIOc+DAAe677z6io6OpUqUK//jHP1i3bp2Ch4hIJaOeDylzhmEwa9Ys/va3v5GZmckf/vAHoqOjad26tbNLExERJ7C75+PSpUsMGzaM119/nVGjRpGRkXFdG6vVynPPPYevry/Nmzdnz549pVKsuJ7ff/+drl27Mnz4cDIzM+ncuTMHDhxQ8BARqcTsDh8RERG0b9+eCRMmEBwczNixY69rM2XKFJ5++mm2bt1KvXr16NKlC5cuXSqVgsV17N27l+DgYFasWIGHhwfTp09n1apV1KqlrdFFRCozk2EYRlEbJyQk0KhRIy5evIjZbObcuXPUr1+fM2fO4OPjk9duy5YttGvXDrgywdDPz49///vftGzZ8qafkZycjMViwWq14uvrW4wfSZzNMAxmzJjB6NGjycrKomHDhixdurRI//+LiIhrsuf7266ej23btuHn54fZfOVE0Tp16uDp6UlsbGy+drnBA8DX1xdfX19uu+22Aq+ZkZFBcnJyvoe4rqSkJJ566ilefPFFsrKyeOaZZ9i/f7+Ch4iI5LErfMTHx1/XZe7j40NCQkKh7zly5Aht27YlIKDg8zgmTZqExWLJe9SrV8+ekqQc2bNnD82bN2f16tV4enry4Ycfsnz5cmrUqOHs0kREpByxK3yYTKa8Xo9cmZmZeHh4FPqejz76iPfee6/QPx87dixWqzXvcfLkSXtKknLAZrPx3nvv8eCDD3LixAkaNWrE7t27GTZsGCaTzmIREZH87FpqGxgYiNVqzfdaamoqgYGBBbb/6quvaNeuHQ0bNiz0ml5eXnh5edlThpQj58+fp0+fPqxfvx6AsLAw5s6dq/k6IiJSKLt6Ptq2bcupU6fIzMwEyBtuCQkJua7tjz/+yIkTJ+jSpUsplCnl0bfffktQUBDr16/Hy8uLOXPmEBUVpeAhIiI3ZFf4CAwMJDQ0lO3btwOwadMmhg4ditlsZurUqRw+fBiAY8eOMW/ePDp06MCvv/7K999/z+zZs0u/enEKm83Gu+++yyOPPEJ8fDxNmjQhNjaWQYMG5RtmybEZ7D52gdVx8ew+doEcW5EXVomISAVm9w6nkZGRjBkzhpiYGJKSkpg8eTIAUVFRNGjQAF9f37wekhkzZuS9L/csD3FtZ8+epVevXmzatAmAXr16MWvWLKpXr56v3caDiYxfe4hEa3reawEWM+M6NSW0WcGTj0VEpHKwa58PR9A+H+XXtm3b6NGjB4mJiXh7e/PRRx/Rp0+f6yaVbjyYSMSi/Vz7i5XbanZ4sAKIiEgFU2b7fEjllJOTw/jx43n00UdJTEykadOm7N27l759+14XPHJsBuPXHroueAB5r41fe0hDMCIilZjCh9zQ6dOn+ctf/sJbb72FzWajX79+7N27l7vuuqvA9rHHk/INtVzLABKt6cQeTyqjikVEpLzTqbZSqM2bN9OzZ0/Onj1LtWrVmD17Nr169brhe86mFB48itNOREQqHvV8yHWys7N5/fXX+ctf/sLZs2e5++672bdv302DB0BdH/NN29jTTkREKh6FD8knPj6eRx99lIkTJ2IYBoMHDyYmJoY777yzSO8PaViLAIuZwvY1NXFl1UtIQ51sKyJSWSl8SJ4NGzYQFBTEjh078PHxISoqisjISLy9vYt8DXc3E+M6NQW4LoDkPh/XqSnubtp2XUSkslL4ELKysnj11Vd5/PHHOX/+PM2bN+e7776jW7duxbpeaLMAZocH42/JP7TibzFrma2IiGjCaWV34sQJunfvzq5duwAYPnw477333nUHCNortFkAHZr6E3s8ibMp6dT1uTLUoh4PERFR+KjE1q5dS+/evbl48SIWi4VPPvmEZ555Jl+bHJtR7ADh7maidaPaZVG6iIi4MIWPSigzM5OxY8cybdo0AFq2bEl0dDS33357vnbaIl1ERMqC5nxUMsePH+fBBx/MCx6jRo1i586dBQaPiEX7r9sw7LQ1nYhF+9l4MNFhNYuISMWi8FGJrFy5kubNmxMbG0vNmjVZvXo106dPx9PTM187bZEuIiJlSeGjEsjIyGDEiBE888wzWK1WWrVqxYEDB+jcuXOB7bVFuoiIlCWFjwru8JFfuPe+ED788EMAXn7lFXbs2EH9+vULfY+2SBcRkbKk8FGBvTZ1DnfdE8ThH3/AzduXus+OY2fNx/jm8Pkbvk9bpIuISFlS+KiALl++zBNhvZn0yhByMtLwuq0pAX0+wLtRyyJNGNUW6SIiUpYUPiqYw4cP06pVK75ctgAw4ds6jFu6T6KKrx9QtAmj2iJdRETKksJHBbJ48WLuu+8+fvjhB9yqWqjb9W1qPtQLk5t7vnZFmTCqLdJFRKSsaJOxCiAtLY0XXniBTz75BIC7Wz5AUsgQqlS/8bDIzSaMaot0EREpCwofLu7QoUN07dqVH3/8EZPJxJtvvkn7HkMJ/3TvTd9blAmj2iJdRERKm8KHC/v8888ZNmwYaWlp+Pv7s3jxYtq1a0eOzSDAYua0Nb3AjcJMXBk+0YRRERFxBs35cEGpqan07t2bvn37kpaWRocOHYiLi6Ndu3aAJoyKiEj5pvDhYv7zn//QsmVLFixYgJubGxMmTGDjxo3ccsst+dppwqiIiJRXGnZxEYZh8PHHH/PCCy+Qnp5OYGAgUVFRPPTQQ4W+RxNGRUSkPFL4cAEpKSkMHjyYqKgoADp27Mj8+fOpU6fOTd+rCaMiIlLeaNilnDtw4ADBwcFERUXh7u7OlClTWLduXZGCh4iISHmkno9yyjAMZs+ezYsvvkhmZib16tUjOjqaBx54wNmliYiIlIjCRzlktVoZMGAAK1asAKBz58589tln1KqlpbEiIuL6FD7KmX379tG1a1eOHz+Oh4cHkyZPpnXn5/n2xGXqXrygCaMiIuLyFD7KCcMw+OCDD3jllVfIysqiQYMGjHr3I6J/9WLmxzF57QIsZsZ1aqqlsiIi4rI04bQcuHjxIk8//TSjRo0iKyuLp59+mvcWb+D97w0SrfnPXzltTSdi0X42Hkx0UrUiIiIlo/DhZHv27KF58+asWrUKT09PZs6cydJly/nn9vgCt0bPfW382kPk2ApqISIiUr4pfDiJzWZj6tSpPPjgg/z22280atSI3bt3M3z4cPb+evG6Ho+rGUCiNZ3Y40mOK1hERKSUaM6HE1y4cIHevXuzfv16AMLCwpg7dy6+vr7AzY+6z1XUdiIiIuWJwocD5NiMvC3OEw/HMfHlIZw6dQovLy9mzJjBoEGDMJn+t4KlKEfd29NORESkPFH4KGMbDyYyfu0hEn5PI3nPCn7/dhEYNm5tcDvrV63k3nvvve49IQ1rEWAxc9qaXuC8DxNXDogLaah9P0RExPVozkcZ2ngwkYhF+zmVcJqzy8bx+44FYNiodtcjuD81hUT3ugW+z93NxLhOTYErQeNquc/HdWqq/T5ERMQlKXyUkRybwfi1h7h84gcSP3+B9F8PYKriRe2OL1D7ib/h5ul9wxUroc0CmB0ejL8l/9CKv8XM7PBg7fMhIiIuS8MuZWT30XP89OVnWP8dBYYNj9p/wK/Lq3jWqQ/kX7FS2Kmzoc0C6NDUP2++SF0fs3Y4FRERl6fwUQZOnz7N0F7PYo39NwDV7m5PrfZDcPO8foLozVasuLuZCg0nIiIirkjho5Rt3ryZ8PBwzpw5g8nDTK2/DKV6s3aFtteKFRERqWwUPorg6qWyhQ19ZGdnM378eCZOnIhhGNx99924PfoiVq+6WrEiIiJyFYWPm8hdKnv1jqPXHu4WHx9Pjx492LFjBwCDBg3i/fffZ/ux34lYtB8T5AsgWrEiIiKVmckwjHJ1QEhycjIWiwWr1Zq346ez5C6VvfYG5caF2eHBcOp7evXqxfnz56levTrz5s2jW7du+a5xs/AiIiLi6uz5/q404aMoQyfXtm8zZUvhZ6zkZJO1N5qE7dEABAUFsWzZMho3blzizxYREXE19nx/V4phl+L0PsQeTyo0eGQnn+P8mn+QEf8TAMOGDWPq1KmYzQVPHtWKFRERkf+p8JuM5Q6dXBskTlvTiVi0n40HEwt8X2FLYNOOxpD42QtkxP+EybMqo/8RyYcfflho8BAREZH8KnT4yN1ltKBxpdzXCttl9NolsEZOFklbPubcv97Blp6Cp39jAvp+wF+ffrb0CxcREanAKvSwy42GTuDGu4xefbhblvUM51ZPITPxCAA+LbpQq20fAmr5aKmsiIiInSp0z8fNdg+9Ubvcw90uHdlF4mcvkJl4BDevatR5+nVqPzoQk7uHlsqKiIgUQ4Xu+Sjq7qEFtcvIyODLuZM498VMADwDm1Cn86tUsdTFX0tlRUREiq1Ch4+rh07s2WX02LFjhIWF8d133wHw0ssv06X/SySl52iprIiISAlV6PCRO3Rizy6jy5cvZ8CAASQnJ1O7dm3mz5/PE0884ciyRUREKrQKPecDrhxLPzs8GH9L/qEVf4uZ2eHBeUMn6enpDB06lK5du5KcnEybNm2Ii4tT8BARESllFbrnI1doswA6NPUvdJfRI0eO0LVrV77//nsAxo4dy9tvv02VKpXi9oiIiDhUpfl2LWyX0SVLljB48GBSU1OpU6cOCxcu5LHHHnNChSIiIpWD3eHj0qVLjB49mpo1a5KamsqUKVPw8vK6rt2aNWv45ptvyMjI4Nlnn6V9+/alUnBpSUtL44UXXuCTTz4BoG3btixevJjAwEAnVyYiIlKx2R0+IiIieOqpp3jqqadYsGABY8eOZdq0afna/Pzzz0yYMIGYmBgMw6BFixasXbuWW2+9tdQKL4mffvqJrl27cvDgQUwmE2+++SZvvPEG7u7uzi5NRESkwrNrwmlCQgLLly+nY8eOAHTs2JHIyEhSUlLytXv//fcJDQ3FZDLh5uZG69atmT17dulVXQLz58+nRYsWHDx4EH9/fzZv3sxbb72l4CEiIuIgdoWPbdu24efnl3eIWp06dfD09CQ2NjZfuy1btlC/fv28540bN2b79u0FXjMjI4Pk5OR8j7Jw6dIlevfuTZ8+fUhLS6N9+/bExcXRrl27Mvk8ERERKZhd4SM+Pp5atfJvyOXj40NCQsIN2xXUJtekSZOwWCx5j3r16tlTUpHNmTOHBQsW4ObmxoQJE9i4cSO33HJLmXyWiIiIFM6u8GEyma47Oj4zMxMPD48btiuoTa6xY8ditVrzHidPnrSnpCJ74YUXCAsLY+vWrfz973/XMIuIiIiT2DXhNDAwEKvVmu+11NTU61aIXNsuJSWl0FUkXl5eBa6WKW1VqlQhOjq6zD9HREREbsyuno+2bdty6tQpMjMzAfKGUkJCQvK1e/TRRzly5Eje86NHj/LII4+UtFYRERGpAOwKH4GBgYSGhuZNHt20aRNDhw7FbDYzdepUDh8+DMCQIUPYvHkzANnZ2cTGxjJw4MBSLl1ERERckd37fERGRjJmzBhiYmJISkpi8uTJAERFRdGgQQOaNGnCvffeS9++fXn55ZfJzMxk+vTp+Pv7l3rxIiIi4npMhmEUdNq80yQnJ2OxWLBarfj6+jq7HBERESkCe76/K/yptiIiIlK+KHyIiIiIQyl8iIiIiEMpfIiIiIhDKXyIiIiIQyl8iIiIiEMpfIiIiIhDKXyIiIiIQyl8iIiIiEPZvb16WcvdcDU5OdnJlYiIiEhR5X5vF2Xj9HIXPlJSUgCoV6+ekysRERERe6WkpGCxWG7Yptyd7WKz2UhISMDHxweTyVSq105OTqZevXqcPHlS58aUId1nx9B9dgzdZ8fQfXacsrrXhmGQkpJCYGAgbm43ntVR7no+3NzcuO2228r0M3x9ffXL7QC6z46h++wYus+OofvsOGVxr2/W45FLE05FRETEoRQ+RERExKEqVfjw8vJi3LhxeHl5ObuUCk332TF0nx1D99kxdJ8dpzzc63I34VREREQqtkrV8yEiIiLOp/AhIiIiDqXwISIiIg6l8CEiIiIOVe42GSsNly5dYvTo0dSsWZPU1FSmTJlS4KzeNWvW8M0335CRkcGzzz5L+/btnVCt6yrKfbZarQwYMICvvvqKRo0aMXv2bFq1auWkil1TUX+fc0VHRxMZGcm2bdscV2QFYM99TktLY86cOdSuXZs77rhDv9N2KOp9fvfdd7HZbLi5uZGWlsY777xT6rteVwYbNmxg3LhxLFu2jAYNGhTYxinfhUYF1KtXL2PlypWGYRjG/PnzjRdffPG6Nj/99JPRsmVLw2azGTk5OUbz5s2NU6dOObpUl1aU+zx27FhjyZIlxr59+4xOnToZdevWNVJTUx1dqksryn3OFR8fbzRr1sx4+OGHHVRdxVHU+3zhwgXj8ccfN/7v//7PkeVVGEW5z2vXrjUGDBiQ97xv375GdHS0w2qsKM6cOWOsWbPGAIzjx48X2MZZ34UVLnzEx8cbZrPZuHz5smEYhnH27FnD29vbSE5Oztdu8ODBxhtvvJH3fOjQocbf//53h9bqyop6n7/55pu8f7ZarYaHh4cRGxvr0FpdWVHvc64RI0YYH330kcKHney5zx06dDC2bt3q4AorhqLe53/84x9GWFhY3vOhQ4cac+bMcWitFUVOTs4Nw4ezvgsr3JyPbdu24efnh9lsBqBOnTp4enoSGxubr92WLVuoX79+3vPGjRuzfft2h9bqyop6n9u1a5f3z7nnCJT12T0VSVHvM8C8efMIDw+natWqji7T5RX1Pq9bt45ffvmF2NhYHn/8ccaOHUtWVpYzSnZJRb3PnTt3Zs2aNaxYsYKLFy9y/vx5evXq5YySXd7NDnhz1ndhhZvzER8fT61atfK95uPjQ0JCwg3bFdRGClfU+3y1I0eO0LZtWwICAsq6vAqjqPf56NGjJCcnExISwqFDhxxZYoVQ1Pu8ePFi7r//fkaNGkX//v1p3rw5hmEwefJkR5brsop6n5s0acKSJUvo0aMHDz/8MP/617/w9vZ2ZKmVhrO+Cytcz4fJZMpL1bkyMzPx8PC4YbuC2kjhinqfr/bRRx/x3nvvlXVpFUpR7nNOTg6zZs1i1KhRDq6u4ijq7/OPP/5ImzZt8PT0pHbt2gwcOJAFCxY4slSXZs9/N9LT01mzZg3Hjh2je/fuZGdnO6rMSsVZ34UVLnwEBgZitVrzvZaamkpgYOAN26WkpFzXRgpX1Puc66uvvqJdu3Y0bNjQEeVVGEW5z7t27SIyMpLatWtTo0YNhg4dys6dO6lRo4aDq3VdRf19zs7OJicnJ+/5PffcQ1JSkkNqrAiKep+3bdvGvn37aN++Pdu2bePAgQNMmzbNkaVWGs76Lqxw4aNt27acOnWKzMxMgLzuo5CQkHztHn30UY4cOZL3/OjRozzyyCOOK9TFFfU+w5W/LZ44cYIuXbo4tMaKoCj3uWXLlhw6dIi4uDji4uJ4++23adGiBXFxcc4o2SUV9ff5nnvu4Zdffsl7XqVKFe68807HFeriinqfV6xYQePGjYErX44TJkxgx44dji22knDWd2GFCx+BgYGEhobmTZjZtGkTQ4cOxWw2M3XqVA4fPgzAkCFD2Lx5M3DlbzOxsbEMHDjQaXW7mqLe52PHjjFv3jw6dOjAr7/+yvfff8/s2bOdWbpLKcp9NpvNNGjQIO+RO6GvsDX9cr2i/j6PHDmS9evXk5GRAVzpdRoxYoTT6nY1Rb3PQUFBHDhwIO99JpOpwL/YyM0Z/z071rjqDNny8F1Y4SacAkRGRjJmzBhiYmJISkrKmwwWFRVFgwYNaNKkCffeey99+/bl5ZdfJjMzk+nTp+Pv7+/kyl3Lze6zr69v3t90ZsyYkfe+6OhoZ5Xskory+ywlV5T73Lp1a9555x1GjBhBkyZNMJvN9O/f38mVu5ai3Od+/fpx4sQJJk+ejJ+fHydPnuS1115zcuWuJzU1lYULFwIwf/58hg8fjp+fX7n4LjQZV8chERERkTJW4YZdREREpHxT+BARERGHUvgQERERh1L4EBEREYdS+BARERGHUvgQERERh1L4EBEREYdS+BARERGHUvgQERERh1L4EBEREYdS+BARERGHUvgQERERh/p/90yF7qyMRqsAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "x1 = np.linspace(0,1,1000)\n",
    "plt.plot(pc,lognorm.cdf(d,params[0],scale=params[2]),'o')\n",
    "plt.plot(x1,x1,color='k')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 63,
   "id": "17bd70c2-a343-4a68-9aa3-568a6f46449e",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "lilfor = np.zeros(len(d))\n",
    "\n",
    "for i in range(len(d)):\n",
    "    lilfor[i] = abs(pc[i]-lognorm.cdf(d[i],params[0],scale=params[2]))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 64,
   "id": "5d71d795-bd7f-4d8f-926e-e8b4248fb2e2",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.07744635678031964"
      ]
     },
     "execution_count": 64,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "max(lilfor)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "324e4884-47b3-406f-a9b5-b0ce3a66176c",
   "metadata": {},
   "outputs": [],
   "source": [
    "w = np.zeros(30)\n",
    "\n",
    "for i in range(len(w)):\n",
    "    uplift = np.loadtxt('./Results/UpliftI_GM'+str(i+1)+'.out')\n",
    "    w[i] = max(abs(uplift[:,1]))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "5c2f7af0-da93-46d3-87a7-e0aac68f4d43",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "w = np.sort(w)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "id": "40f0304f-0b0e-4ec3-8408-e7437d5a75df",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.07753791999999998"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "np.mean(w)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "id": "f4f1bb3f-56f2-46a9-b4c5-d69f043df917",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.07929335"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "np.median(w)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "id": "100bd885-08f6-4ca4-bc25-65f580f02b72",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "pcw = np.linspace(0,1,len(w))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "id": "b5d432c6-5fe7-4cd4-97c5-6605e9d3342e",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "paramsw = lognorm.fit(w,floc=0.0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "id": "dbdcdfb1-de5c-45f6-9320-0bc3c2ed9a28",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(0.4893180373796253, 0.0, 0.06947457022237483)"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "paramsw"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "id": "e1d725d8-82e8-4bfc-8ce0-9461c657e076",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x2e5c31550d0>]"
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiUAAAGdCAYAAADNHANuAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABATElEQVR4nO3deXhU5f3+8fdMtgFNwhIwCSBBBBQRJWoU6wIS/AYVEAtaFbGKqKAo2tYaq1KEClR+Lq0KFqtVUSiLrSAVqWCwdSEioAISVtkS1uBkgWSSmfP74zSBIZNkJplJZrlf15XLWZ458zyFntyc5fOxGIZhICIiItLMrM09ARERERFQKBEREZEgoVAiIiIiQUGhRERERIKCQomIiIgEBYUSERERCQoKJSIiIhIUFEpEREQkKEQ39wRO5XK5yM/PJz4+HovF0tzTERERES8YhkFxcTGpqalYrQ075hF0oSQ/P59OnTo19zRERESkAfbs2UPHjh0b9NmgCyXx8fGAuaiEhIRmno2IiIh4o6ioiE6dOlX/Hm+IoAslVadsEhISFEpERERCTGMuvdCFriIiIhIUFEpEREQkKCiUiIiISFBQKBEREZGgoFAiIiIiQUGhRERERIKCQomIiIgEBYUSERERCQpBVzxNREREAsPpMsjdWcjB4jLax9vI6NKGKGvw9JlrUCj56KOPmDhxIvPnzyctLc3jmMWLF7NixQrKy8sZPnw4mZmZjZmniIiINMKyDQVMWrKJAntZ9WspiTYmDu5JVq+UZpzZCT6HkoMHD1JZWcnXX39d65jNmzczZcoUVq9ejWEYXHzxxSxZsoQOHTo0arIiIiLiu2UbChg7Zy3GKa/vt5cxds5aZo5MD4pg4vM1Je3bt+f666+vc8yLL75IVlYWFosFq9VK3759mTlzZoMnKSIiIg3jdBlMWrKpRiABql+btGQTTpenEU2rQRe6Wq11f2zlypV07ty5+nm3bt1YtWqVx7Hl5eUUFRW5/YiIiIh/5O4sdDtlcyoDKLCXkbuzsOkmVYuA3H2zb98+2rRpU/08Pj6e/Px8j2OnTp1KYmJi9U+nTp0CMSUREZGIdLC49kDSkHGBFJBQYrFYsNls1c8dDgcxMTEex2ZnZ2O326t/9uzZE4gpiYiIRKT28bb6B/kwLpACcktwamoqdru9+nlxcTGpqakex8bFxREXFxeIaYiIiES8jC5tSEm0sd9e5vG6EguQnGjeHtzcAnKkZMCAAWzZsqX6+bZt2+jfv38gvkpERETqEGW1MHFwT8AMICerej5xcM+gqFfSoFBiGIbbfwFmzJhBXl4eAPfffz+ffPIJAJWVleTm5jJmzJjGzlVEREQaIKtXCjNHppOc6H6KJjnRFjS3A0MDTt+UlJTwzjvvAPDWW2/x4IMPkpSUxNy5c0lLS6NHjx5ccMEF3HXXXfz617/G4XDwwgsvkJyc7PfJi4iIiHeyeqUwsGdyUFd0tRgnH+4IAkVFRSQmJmK320lISGju6YiIiIgX/PH7W71vREREwlCw97nxRKFEREQkzIRCnxtPAnL3jYiIiDSPqj43p1Zxrepzs2xDQTPNrH4KJSIiImEilPrceKJQIiIiEiZCqc+NJwolIiIiYSKU+tx4olAiIiISJkKpz40nCiUiIiJhoqrPTW03/low78IJhj43niiUiIiIhIlQ6nPjiUKJiIhIGAmVPjeeqHiaiIhImAmFPjeeKJSIiIiEIsOAkgPgKIW2XWu8HWW10Ldr22aYWMMplIiIiPhBQHvNuFxwdCcUfAv7v4OC78z/lh6CbtfC7Qv88z3NTKFERESkkfzaa6bSAYd+OBE89n8P+zeAo7jmWIsVKoOz5khDKJSIiIg0QlWvmVMLt1f1mqnz4lLDgMIdsO+bEz8F34GzvObYaBuccR4knw/JvSHlAmjfE2Jb+n1NzUWhREREpIHq6zVjwew1M7Bnsnkqp+Qg7FvrHkLKfqr5YVviieCR3BtSekPbbhAV3r+2w3t1IiIiAVRXrxkLLrpb9nJJSR6Fb/+NdkfXg313zYFRcWbo6HARdLgYOqRDm7PAEtx3ygSCQomIiEgDndxDJg4HF1i2c7E1j0useVxk3UqC5Zj55o9VoyzQrsf/Aki6+d/250F0bFNPPSgplIiIiDREeQndi1bzePQiMqyb6WXZSazF6Tak1IhjrasbaX2uodMF/c0QYktspgkHP4USERERb1Qchz2rYed/YOdnkL+Wc12VnHvSb9IDRiu+dvVgjasHa1zd2Wx0pl3iafz3xmsgyAuXBQOFEhEREU8qHbBvzYkQsjcXnA73Ma3OZG+rS3hhSztyjXPYY7SjqstMKPSaCTYKJSIiIlWO/gjbVpg/O1eBo8T9/fhU6HIldLkK0q6E1p3pCAzcUMAXSzbBSRe9Jje0TkkEUygREZHI5TgGuz6HbZ+YP0e2ub/fMumkEHKVWc7dw10xodprJtgolIiIiF8EtMy6Px3ZDnkfmSFk1xfuhcosUdDpUjh7AJydadYIsVq92mwo9poJNgolIiLSaH4ts+5vLifsXQN5/zLDyOE89/cTO50IIV2u0t0xzUihREREGqVRZdYDxXEMduRA3lLY8rHZuK6KNRo6/wy6Z5lBJKlbRBYqC0YKJSIi0mA+l1kPpLIi80jIpg9g+wr3RnVxidBtIPQYZP5XR0OCkkKJiIg0WF1l1sEMJgX2MnJ3FgbmeovqIPJP8xqRk2/ZTTwTzrnODCJnXq6qqSFAoURERBrs5DLr/hjnlTI75C2Djf8wj4icHESSukPPG6HnEDijl07LhBiFEhERabD28Ta/jqtVRRls/Ri+mw9bl3sOIucNg/bnKoiEMIUSERFpsIwubUhJtLHfXubxuhILZhGxjC5tfN+4ywW7v4Dv/g4bP4By+4n3krqbIaTnjQoiYUShREREGizKamHi4J6MnbMWC7gFkwaXWT+42Qwi3y8A+54Tryd0gPNHQO+boX1PBZEwpFAiIiKNktUrhZkj02vUKfGpzHqZHb5fCOvegfx1J16PS4CeQ6H3LeZtvF4WMpPQpFAiIiKN1qAy64ZhVlRd9w5s/CdUHjdft0ZDt2vNIyLdsyCmRZOsQZqfQomIiPiF12XWiw/At+/BujnuvWbanQPpo8yjIqclBW6iErQUSkREIlCT96kxDNj5GXw9Gzb/Cwyn+XrMadDrJki/EzperOtEIpxCiYhIhGnSPjVlRfDtPPj6dfeeMx0zzKMi5w2DuNP9+50SshRKREQiSJP1qTn4A+TONu+icZSYr8WcBhf8Ai65B87o2fjvkLCjUCIiEiEC3qfG5TRLvn81E3b998TrSd3hkjFmILElNHD2EgkUSkREIkTA+tQ4jsH6d+GrV6Fwh/maxQo9roOMe6HLVbpWRLyiUCIiEiH83qem+ADk/gXW/BWOHzVfsyXCxXebp2gSOzZwphKpFEpERCKE3/rUHPwBvnzZ7ENT1YOmVWfo+wBceLsuXJUGUygREYkQje5Ts3cNfDYDtnx04rWOl0DfB+HcwWCNCsS0JYIolIiIRIgG9akxDNj1OXz2HOzIOTH63Bug73g489ImmbtEBoUSEZEI4nWfGsOAbZ+YR0b2fGW+Zokyq61e8Qi0694Ms5dwp1AiIhJh6uxTYxiwZRnkTIOC9eYHomKhz0j42cPQOq05py5hTqFERCSE+Ks8fI0+NYYB21fCyimw7xvztZiWcNFdcPl4SPBzpVcRDxRKRERCRMDKw//4OXz6B/PaEYDoFnDpvXD5Q2qMJ01KoUREJAQEpDz83jXmkZEdn5rPo+LMGiNXPALxZ/hl3iK+UCgREQlyfi8Pf3grfPJ72Pyh+dwabTbHu/LXkNjBfxMX8ZFCiYhIkPNbefiSQ7BqGqx5EwynWQr+gtvg6t/oAlYJCgolIiJBrtHl4R2l8OWr8PmLJzr2ds+CzEnQ/hz/TFLEDxRKRESCXIPLw7ucZqO8T5+F4gLztdQ+MHAydLnSz7MUaTyFEhGRINeg8vA/fg4f/RYOfG8+b3UmDJgI590EVmtTTFvEZ/qbKSIS5KrKw8OJcvBVapSHt++FhXfD364zA4ktEa79Azy4Bs4frkAiQU1/O0VEQkBVefjkRPdTNMmJNvN24B6tYdVz8PIlsGERYDELn41fB5c/CNFxzTNxER/o9I2ISIjwWB4+rTVRW5bCK7+Dn3aZA8/sC4OmQ8oFzTthER/5HEpKS0t57LHHaN26NSUlJUyfPp24uJoJ/Nlnn8XlcmG1Wjl27BiTJ0/GYvG9FLKIiJzgVh6+cCfMvddsnAcQnwrXToZePwftbyUE+Xz6ZuzYsWRmZjJlyhTS09PJzs6uMebDDz9k586dPPnkkzzxxBPk5+czf/58v0xYRCRSOF0GX24/wgfr9/Hl9iM4Xf+7zNVZAf95Hl69zAwkUbFw5a/gwa/N60YUSCRE+XSkJD8/nwULFvCXv/wFgEGDBnH//fczadIk4uPjq8f98MMPFBcXVz9v0aIFdrvdT1MWEQl/tfW5eaFvGZdtnAyHfjBfTLsSbngBkro100xF/MenIyU5OTkkJSVhs5kXWrVr147Y2Fhyc3Pdxg0ZMoTFixezcOFCjh49yuHDh7njjjs8brO8vJyioiK3HxGRSFbV5+bkQJJICQ+X/pnLcm4zA0nLtnDjLLhziQKJhA2fjpTs27ePNm3auL0WHx9Pfn6+22s9evTgvffe47bbbuPqq69m0aJFtGjRwuM2p06dyqRJk3yctohIeKrZ58bgButX/D7mLZIs5j/aFlsHcP242USdXkdJeZEQ5NOREovFUn2UpIrD4SAmJqbG2LKyMhYvXsz27du59dZbqays9LjN7Oxs7HZ79c+ePXt8mZKISFg5uc9NO35iVsyLvBz7Z5IsRWx1dWBE+dM8dGw0uQeaeaIiAeBTKElNTa1xbUhJSQmpqalur+Xk5LBmzRoyMzPJyclh3bp1PP/88x63GRcXR0JCgtuPiEikMvvXGAy1/pflcY+RFfU1FUYUL1bexHWOqXxtnHPSOJHw4lMo6devH3v37sXhcABUn7bJyMhwG7dw4UK6dTPPcaampjJlyhQ+++wzf8xXRCSsdYgqYnbM87wU+yqtLSVscKUxxDGFFyuHU3HSGXdv++GIhBKfj5RkZWWxatUqAJYvX864ceOw2WzMmDGDvLw8AC688ELWrVtX/TmLxVIjuIiIyEkMA76dx0VL/4+BUd/gMKKYUTGCGx3P8IPRuXqYBfMuHLc+NyJhwufiabNmzeLxxx9n9erVFBYWMm3aNADmzp1LWloaPXr04O6772b37t1MmzaNpKQk9uzZwxNPPOH3yYuIhIXSw7D4IchbigWwtz6PW/bfQZ5xplsDvhp9bkTCjMUwDE9NJ5tNUVERiYmJ2O12XV8iIuFv67/hn+Og9KBZBK3f43D5wyz74ZDHOiUTB/ckq1dKM05YxDN//P5W7xsRkebgOAb/fhq+nm0+b3cu/Hw2JJ8P1NLnpksbHSGRsKZQIiLS1PLXw/v3wmHzOjwuHQuZv4cY94tX3frciEQAhRIRES84XUbjj1q4nPDFn2DlH8BVAacnw42vwtkDAjNpkRCjUCIiUo/a+tD4dH1H8QF4/x7Y+b/yCOcOhsF/gpa6i0akis9dgkVEIomnPjQA++1ljJ2zlmUbCurfyI4cmHWFGUhiToMhL8PN7yiQiJxCoUREpBY1+9CcUPXapCWbcLpquYnR5YRPn4W3bzTvrmnfE+7NgfQ7wKILVkVOpVAiIlKLk/vQeGIABfYycncW1nyzqADeHgqrppsj00fBPSugXfeAzVck1OmaEhGRWnjbX6bGuG0rzLtrjh2G2NPhhheh9wj/T1AkzCiUiIjUwtv+MtXjXC7zyEjV0ZEzesGIv0FSt4DNUSScKJSIiNQio0sbUhJt7LeXebyuxAIkV/WhOX7UPDqydbn55kW/hKxpENOiCWcsEtp0TYmISC2irBYmDu4JnOg7U8WtD82hTfCX/mYgibbBjTNh8EsKJCI+UigREalDVq8UZo5MJznR/VROcqKNmSPTyTI+h9cz4ehOSDwT7v4YLrytmWYrEtp0+kZEpB4e+9B0TiBqxe/hy5fNQWf1h+FvqPaISCMolIiI1OHU8vI39E4lquwovHvTieqsVzwC1zwF1qjmnaxIiFMoERGphafy8pfFH+avcf+P00p2mbf73vgq9BzajLMUCR8KJSIiHlSVlz/5rpurrN/ysuPPnFZxjOMtU2lx50I447xmm6NIuNGFriIip6hZXt7grqiPeDPmjyRYjvG1qzs/r5yCs13PZpylSPhRKBEROcXJ5eWjqeTZ6L8yMeYdoiwGCyqv4nbH79hUZPNcXl5EGkynb0RETlFVNj6BEmbFvMjlUZtwGRamVt7KbOf1VFUp8bYMvYh4R6FEROQU7eNtdOAQb8b+ke7WfZQYNh6qeJCVrvQa40TEfxRKREROkWHbzQe2iSTxE/uN1tzt+A2bjLTq993Ky4uI3yiUiIicbOu/iZp/J0mUstnVibsdj5FP2+q33crLW08tPi8ijaELXUVEqqx5E967BSpK4ax+7B32PkZiB7ch1eXle6U00yRFwpeOlIiIGAasnAz/+X/m8wtug8EvkRkdS/8LurmXl+/SRkdIRAJEoUREIpuzEj6cAOveMZ9f/Tj0exwsZvCIslro27Vt7Z8XEb9RKBGRyFVRBotGw+YPwWKFwS/hvPAOcnfoyIhIc1AoEZHIVFYE826DH/8DUXEw/A2WOS9i0vSVbr1uUhJtTBzcU9eQiDQBXegqIpGn5CD87XozkMTGw8hFLHNexNg5a90CCcB+exlj56xl2YaCZpqsSORQKBGRyHJ0F7zxf7D/O2iZBL/8EGfnK07pdXNC1WuTlmzC6fI0QkT8RaFERCLHwc1mICncAa3OhNHLIfVCt143nhhAgb1MvW5EAkzXlIhIZCj4Dt65EY4dgXbnwh3/gATzOhFve9io141IYCmUiEj42/cNvDMMyuyQcqEZSFqeKBHvbQ8b9boRCSydvhGR8LbrS3hrqBlIOmbAnYvdAglARpc2pCTaqO3GXwvmXTjqdSMSWAolIhK+dqyCOTeBoxjSrjSPkNgSawyLslqYOLgnQI1gol43Ik1HoUREwtOW5fDuCKg4Bl0HwG3zIe70Wodn9Uph5sh0khPdT9Go141I09E1JSISfvI+gr/fAa4K6HE9jHgTouPq/VhWrxQG9kxWrxuRZqJQIiLhJW/ZiUBy3jC4aTZExdT7MafLcAsjN/ROVRgRaWIKJSISPrYsh/knB5LXIar+3dyyDQVMWrJJ5eVFmpmuKRGR8LD1E/j77eB0QM+hPgUSlZcXCQ4KJSIS+rZ9YjbXczrg3MHw8796FUicLkPl5UWCiEKJiIS27Sth3u3gLIdzboDhb3p1DQmg8vIiQUahRERC187PYO6tUFkGPa7zKZCAysuLBBuFEhEJTXty4b1fmIGkexaMeAuiY33ahMrLiwQXhRIRCT0F38G7w6GiFM7qDze/7XMgAZWXFwk2CiUiEloObTnRXK/TZfCLd70qjOaJysuLBBeFEhEJHUd3wdtD4dhhSLkAbp8Psac1apMqLy8SPFQ8TURCQ1EBvD0EivOh3Tkw0nNzvYZQeXmR4KBQIiLBr/QIvHMjHP0RWqfBHf+E09r69SuirBb6dvXvNkXENwolIhLcykvMi1oPbYb4VBj1ASTUf0rl1F42OvIhEvwUSkQkeFU6zF42+WuhRRszkLROq/dj6mUjEpp0oauIBCeXCz54wKzYGtMSbl8A7brX+zH1shEJXQolIhKc/v0UfD8frNFmHZKOF9f7EfWyEQltCiUiEnw+/xN8+bL5eMjL0G2gVx9TLxuR0KZQIiLB5dt55lESgIHPwIW3ev1R9bIRCW0KJSISPLZ+Yl5HAtD3Qbj8IZ8+rl42IqFNoUREgsPeb8w7bVyVcP7NMHAyWHy7hVe9bERCm8+3BJeWlvLYY4/RunVrSkpKmD59OnFxnvtOHDt2jNdee422bdvSvXt3LrvsskZPWETC0NEf4b2boeIYdL0Ghr4CVt//zVTVy2bsnLVYwO2CV/WyEQl+Pv+/fuzYsWRmZjJlyhTS09PJzs72OK6wsJARI0Zw4403MmrUKAUSEfHs+FF4d4TZzya5d4M7/lZRLxuR0GUxDMPre+Py8/Pp2rUrR48exWazcejQITp37syBAweIj493G3vttdfyxBNP0K9fP58mVFRURGJiIna7nYSEBJ8+KyIhptIBc26CH/8DCR3gnhVeVWv1hiq6ijQtf/z+9ulISU5ODklJSdhs5r9A2rVrR2xsLLm5uW7jPvzwQ7Zu3Upubi7XXXcd2dnZVFRUeNxmeXk5RUVFbj8iEgEMA5Y8bAaS2NPhtr/7LZDAiV42Qy/sQN+ubRVIREKAT9eU7Nu3jzZt3C8Qi4+PJz8/3+21d999l0svvZQJEyYwevRo+vTpg2EYTJs2rcY2p06dyqRJkxowdREJaZ89B9++B5YoGPEWJJ/v1cd0BEQkfPl0pMRisVQfJanicDiIiYlxe23jxo1cccUVxMbG0rZtW8aMGcPbb7/tcZvZ2dnY7fbqnz179vi4BBEJOd/Nh0//YD6+fgZ0y/TqY8s2FHDF9JXcOvsrHp63nltnf8UV01eqdLxImPAplKSmpmK3291eKykpITU11e21yspKnE5n9fPevXtTWOi5gmJcXBwJCQluPyISxn78/EQtkssfgovv9upj6mkjEv58CiX9+vVj7969OBwOgOrTNhkZGW7jevfuzdatW6ufR0dHc8455zR2riIS6g5vg7/fDk4H9BwKmd6dulVPG5HI4PORkqysLFatWgXA8uXLGTduHDabjRkzZpCXlwfAww8/zNKlSykvLwfgiy++YPz48X6euoiElONHzVokx49Cx0tg2Gte1yJRTxuRyOBz8bRZs2bx+OOPs3r1agoLC6svXp07dy5paWn06NGDvn37MnnyZMaPH0+PHj2w2WyMHj3a75MXkRDhrIQFv4TC7ZB4JvxiLsS08Prj6mkjEhl8DiVJSUm8/vrrNV7/5ptv3J6PHDmSkSNHNnxmIhI+lv8OduRAzGlw61w4vZ1PH1dPG5HIoN43IhJY37wFq2eZj296DZJ7+bwJ9bQRiQwKJSISOLu+gKW/Mh/3fxLOHdygzVT1tAFqBBP1tBEJHwolIhIYR3fB30eCqwLOuwmu+nWjNqeeNiLhz+drSkRE6lVeAvNug2NHIOVCs+uvpfFHMbJ6pTCwZ7IquoqEKYUSEfEvlwv+cR8c2ACntYdfvAexLX3eTG3l5Kt62ohI+FEoERH/ypkKmz+EqFgzkCR28HkTyzYUMGnJJrfaJCmJNiYO7qnTNCJhTNeUiIj/bF4Kn/3RfDz4T9DpEp83oXLyIpFLoURE/OPwVnj/PvPxpWPhwlt93oTKyYtENoUSEWm88mKYdzs4iqHzz+DayQ3ajMrJi0Q2hRIRaRzDgH+Og8N5EJ8CI/4GUTEN2pTKyYtENoUSEWmcz1+EHxaDNQZufgdOb9/gTamcvEhkUygRkYbbvhJWPGM+vu65Bl3YejKVkxeJbAolItIwR3fBwrvBcEGfO+CiXzZ6kyonLxLZFEpExHcVx80S8sePQmo6XDfDLxVbQeXkRSKZiqeJiG8Mw2yyt/87aJkEt7wDMf69xkPl5EUik0KJiPhm3Tuw/l2wWGHEm5DYMSBfo3LyIpFHoUREvLf/e/jXb8zH1zwFXa4K2FfV1vtGRMKXQomIeKesCOaPgsoy6HYt/GxCwL5KvW9EIpMudBWR+hkGLH4QCndAYicY9hpYA7P7UO8bkcilUCIi9Vv9Gmz6wCyQNuJv0DIwdULU+0YksimUiEjd9q6B5U+aj6+dAh0vDthXqfeNSGRTKBGR2h0rhAW/BFcF9BwKl94X0K9T7xuRyKZQIiKeuVzwj/vAvgfanAVD/uy3Amm1Ue8bkcimUCIinn3+ImxdDlFxcPPbYEsM+Feq941IZFMoEZGadq+GlVPMx9c9B8nnN8nXqveNSGRTKBERd8ePwqLRYDjh/BGQPqpJv169b0Qil4qnicgJhgGLHzKvI2ndBa5/PuDXkXii3jcikUmhREROWPMG/LDYrEcy/A2wJfht076WjVfvG5HIo1AiIqYDG+HjJ8zHmROhQ7rfNq2y8SLiDV1TIiLgOAYL7jL72pw9EC57wG+bVtl4EfGWQomIwLLH4XAenH4G3DjTb31tVDZeRHyhUCIS6TYsgrVvARa46S9weju/bVpl40XEFwolIpHs6I+wZIL5+MpH4ax+ft28ysaLiC8USkQilbMCFo6G8iLodCn0y/b7V6hsvIj4QqFEJFJ9+gfYt8YsH//z1yEqxu9fobLxIuILhRKRSLTzP/DfF83Hg/8Erc4MyNeobLyI+EKhRCTSHP8J/nE/YECfkXDejQH9OpWNFxFvqXiaSKRZ+iso2muWkc+a3iRfqbLxIuINhRKRSPLdfNiwECxR5nUkcac32VerbLyI1EehRCRS/LTbPEoCuK56jNXlXTi4fp+OWohI0FAoEYkELie8fx+UF3G0zYUM/jKdvR9/Vf22+tCISDDQha4ikeDzF2H3F1RGt2RowZ3sLapwe1t9aEQkGCiUiIS7fWvh02cBmMrd7DbOqDFEfWhEJBgolIiEM0cpvD8GXJUcOXMQfy3pW+tQ9aERkeamUCISzpY/CUe2QXwKq3s9Rc0SZjWpD42INBeFEpFwlfcRrHnDfHzjTFq3TfbqY+pDIyLNRXffiISjkoPwwYPm474PQtf+ZLgMUhJt7LeX4emqEQtmlVX1oRGR5qIjJSLhxjBg8UNw7DC0Pw8GPA2oD42IBD+FEpFws/5d2PIRRMXCz2dDdFz1W+pDIyLBTKdvRMLJT7vho8fNx/2fgDPOqzFEfWhEJFgplIiEAafLIHfHYc5edg/tHMUYHS/FcvlDtY5XHxoRCUY6fSMS4pZtKOCK6StZ9uZk2h1ezTEjjpsPjmLZpoPNPTUREZ8olIiEsGUbChg7Zy22op08Hj0XgKmVt7KmqLXKxotIyFEoEQlRTpfBpCWbsOLk+ZiZtLA4+I+zF3OcmSobLyIhSaFEJETl7iykwF7GfVFL6GPdRpHRkscq7sP43/+tVTZeREKNQolIiDpYXMa5ll1MiF4EwO8rRlFAzYtXVTZeREKFz6GktLSUBx54gCeffJIJEyZQXl5e5/h58+bRr1+/hs5PRGpxRksLz8e8SqzFycfOi3nfdaXHcSobLyKhwudQMnbsWDIzM5kyZQrp6elkZ2fXOjY/P58//OEPjZqgiHiWsesvnGvdw2EjgScqRnNqnVYLkKKy8SISQnwKJfn5+SxYsIBBgwYBMGjQIGbNmkVxcbHH8dOmTWPs2LGNn6WIuNu9GusXLwHwu4rRFJLo9rbKxotIKPIplOTk5JCUlITNZh4ObteuHbGxseTm5tYYO3v2bEaOHEnLli3r3GZ5eTlFRUVuPyJSB0cp/PN+MFzQ+xcMu+0+lY0XkbDgU0XXffv20aaN+6Hg+Ph48vPz3V7btm0bRUVFZGRksGnTpjq3OXXqVCZNmuTLNEQi278nQuEOiE+FQdPJatFKZeNFJCz4FEosFkv1UZIqDoeDmJiY6udOp5NXX32V5557zqttZmdn8+ijj1Y/LyoqolOnTr5MSyRy7FgFX882Hw99GVq0AlQ2XkTCg0+hJDU1Fbvd7vZaSUkJqamp1c+/+OILZs2axRtvvAGYocXhcNCqVSt++umnGtuMi4sjLi6uxusicoryElj8oPn4ortwnnUNuduP6OiIiIQNn0JJv379uPfee3E4HMTGxlaftsnIyKgec8kll7idslm4cCELFy5k3rx5fpqySIT6ZKLZBTjxTP7d8QGenr6SAvuJGiQpiTYmDu6p60hEJGT5dKFramoqWVlZrFq1CoDly5czbtw4bDYbM2bMIC8vD5vNRlpaWvVP1YWxaWlpgZi/SGTY+Rl8/ToAub0nce/ft7gFEoD99jL1uxGRkOZznZJZs2bx97//nSlTpvDdd99V1yGZO3cu33//vd8nKBLxykvgA/O0jeuiu3g4NxFP3WzU70ZEQp3FMIyg2nsVFRWRmJiI3W4nISGhuacj0vz+9RvI/QskdiI360NufmtjvR+ZO+YyXfgqIk3KH7+/1ftGJJj9+F8zkAAM+RMF5TF1j/8f9bsRkVCkUCISrByl8MED5uP0O6HrNV73sVG/GxEJRQolIsFqxTNw9EdI6AjXTgEgo0sbUhJt1Hbjr/rdiEgoUygRCUY/fg6rZ5mPh7wENvP8bJTVwsTBPYFT2++p342IhD6FEpFg4zh20mmbUXB2ptvbWb1SmDkyXf1uRCTs+FQ8TUSawIpn4OhOSOhQfdrmVFm9UtTvRkTCjkKJSDDZ9cWJ0zaD/wS2xFqHqt+NiIQbnb4RCRaOYxgfPAAY7DrzJr609lERNBGJKDpSIhIkflzwOGmFOygw2jB4yyCKtnylfjYiElF0pEQkCHyVs5Qzt7wNQHbFaIo4DVA/GxGJLAolIs3MWX6M1FW/wmoxWFB5FTmuPtXvqZ+NiEQShRKRZnbgg6c40yhgv9GayZV31HjfAArsZeTuLGz6yYmINCGFEpHmtHs1KZv+CkB2xT3Vp208UT8bEQl3CiUizaXiOHzwABYMFjqv4tOTTtt4on42IhLuFEpEmsunz8KRrRinJ/OXFveon42IRDyFEpHmsPcb+PJlACyDX+TRIRnm41OGqZ+NiEQShRKRplZZDh+MA8MF598MPQapn42ICCqeJtL0Vv0RDm2G09rBoOnVL6ufjYhEOoUSkSbidBls+OY/nP/fF7ACzutmENXS/ToR9bMRkUim0zciTWDZhgL6TfuYmCUPYjWcfOi8lCsWx6tSq4jISRRKRAJs2YYCxs5Zy7DSBfS07qLQOJ2JFb9UCXkRkVMolIgEkNNlMGnJJrpbdvNg9D8A+H3FLzlCokrIi4icQqFEJIBydxZy0F7KczGvEWtxstx5EYtdfavfVwl5EZETFEpEAuhgcRn3Ri2lt3UndqMlv6u4m5rVSFRCXkQEdPeNSEB1du0lK3oRAM9UjOIQrT2OUwl5ERGFEpHAcTm54JvfYbFU8KnzAha5rqwxxIJZIE0l5EVEdPpGJHC+moll39dURp/G7yruwXLKaRuVkBcRcadQIhIIR7bDyskARA96lqdHDlQJeRGReuj0jYi/uVyweDxUlkGXqyH9TrIsFpWQFxGph0KJiL99/Trs+hxiToMhfwaLGTxUQl5EpG46fSPiR84jO3H+eyIAO/r8Bmfimc08IxGR0KFQIuIny77PZ+3Lo4iqPMZq1zkM+Oxsrpi+UmXkRUS8pFAi4gfLNhTw2bwZXGJ8x3Ejlscq7sXAqv42IiI+UCgRaSSny2Dm4s/Ijn4PgBmVN7PLSAZQfxsRER8olIg0Uu6OIzxy/BXiLcf5xtWNN51Zbu+rv42IiHd0941II8VtnEffqG8pN2J4rOJeXLVkffW3ERGpm46UiDRGUQG9N04H4IXKn7Pd6FDrUPW3ERGpm46UiDSUYcDSR4l2FLHJ0pXXndd7HKb+NiIi3tGREpGG2rAI8v4F1hgKM1/ASRSn1mdVfxsREe8plIg0RMlB+NdvzMdXP8YVP7uamSPT1d9GRKQRdPpGxFeGAUsmwPFCSD4frngEgKxeKepvIyLSCAolIvVwugy3oHFpySdY85aCNQZunAVRMdVj1d9GRKThFEpE6rBsQwGTlmyiwG7eztueo3xi+y0JAP1+C8m9mnV+IiLhRKFEpBbLNhQwds5aTtRhNZga8zoJlPCdqwsFrX7B/zXj/EREwo0udBXxwOkymLRkEycXhh8e9RkDotZRbkTz64qx/H7pFpWOFxHxI4USEQ9ydxZWn7IBSOYIT0e/DcALlcPZYnRU6XgRET9TKBHxwL0kvMH0mNkkWI6zznU2s08qkqbS8SIi/qNQIuLBySXhb4nK4eqo7yg3Yvh1xX04ifI4TkREGkcXuop4kNGlDSmJNqLse3gyeg4Az1XeXN3bRqXjRUT8T6FExIMoq4WJN5zL6fOfIt5ynDWu7rzhHASodLyISKAolIjUIqvsI4jaSBmx/KbiPlz/O9uZnGhj4uCeKh0vIuJnCiUinhTuhOVPARD7f5N4tv0wlY4XEQkwhRKRU7mc8MEDUFEKZ16O9dL76WvVNeEiIoGmUCJh59ReNT4f2fjyFdj1OcScBje+AgokIiJNQqFEwsqpvWoAUny5BmT/Blg52Xyc9Sy0OStAMxURkVPpn4ASNqp61ZwcSAD228sYO2ctyzYU1L2BynL4x33gdED3LEi/M4CzFRGRU/kcSkpLS3nggQd48sknmTBhAuXl5TXG2O12RowYQUJCAn369OGrr77yy2RFauOpV02VqtcmLdlUd6+aT5+FAxugZVsY8mew6GJWEZGm5HMoGTt2LJmZmUyZMoX09HSys7NrjJk+fTo33XQTn376KZ06dWLo0KGUlpb6ZcIinpzaq+ZUBtTdq2bXF/D5S+bjwS/B6e39P0kREamTT6EkPz+fBQsWMGiQWURq0KBBzJo1i+LiYrdxmZmZ3HrrrVx00UXMmTOHo0ePsmnTJv/NWuQU3vag8TiurMg8bYMBF94O5w727+RERMQrPoWSnJwckpKSsNnMfh/t2rUjNjaW3Nxct3HXXHNN9eOEhAQSEhLo2LGjx22Wl5dTVFTk9iPiK2970Hgc93E2/LQbEs+ErGl+npmIiHjLp1Cyb98+2rRx7/URHx9Pfn5+rZ/ZsmUL/fr1IyXF850PU6dOJTExsfqnU6dOvkxJBDjRq6a2q0AsmHfh1OhVs3kprJtjjhg2C2wJAZ6piIjUxqdQYrFYqo+SVHE4HMTExNT6mVdeeYXnnnuu1vezs7Ox2+3VP3v27PFlSiLA/3rVDO4JUCOY1NqrpuQQLH7IfHz5eEj7WcDnKSIitfMplKSmpmK3291eKykpITU11eP4jz/+mGuuuYYuXbrUus24uLjqUzxVPyINkdUrhZkj00lOdA/OyYk2Zo5Md69TYhiw+EE4dhjanwfXPNnEsxURkVP5VDytX79+3HvvvTgcDmJjY6tP22RkZNQYu3HjRnbv3s2YMWP8M1MRL2T1SmFgz+T6K7p+/TpsWQZRsXDTaxAd1zwTFhGRaj4fKcnKymLVqlUALF++nHHjxmGz2ZgxYwZ5eXkAbN++ndmzZzNw4EB+/PFHvv32W2bOnOn/2Yt4EGW10LdrW4Ze2IG+XdvWDCQHNsLHvzMfD3wGks9v+kmKiEgNPpeZnzVrFo8//jirV6+msLCQadPMuxXmzp1LWloaCQkJ9OvXj7179/LSSy9Vf27evHn+m7VILerte1NxHBaOBmc5nD0QLr2/+SYrIiJuLIZh1FHisukVFRWRmJiI3W7X9SXiE6/63iz9lXnq5rT2MPYLOL1dM81WRCS8+OP3t3rfSFjwqu/N5n+ZgQRg2EwFEhGRIKNQIiHPm743ryz+L8YHD5hP+j4IZ2c21fRERMRLPl9TIhJs6ut7Ay5+e/wFLFGFkNwbBjzdZHMTERHv6UiJhLz6+t7cF/UhV0RtpDKqBQx/Q7f/iogEKYUSCXl19b1Jt2zh19HzAdh1ydOQ1K2ppiUiIj5SKJGQV1vfm1YU8+fYPxNtcfFv689IG6jbf0VEgplCiYQ8T31vLLj4fzGz6GA5wk5XMgx+iago/XUXEQlm2ktLWDi17829UUsZELWOcmLJv3YWA/votI2ISLDT3TcSNqr63vyw+mN6Lp8PBsTc8Bw/u7h/c09NRES8oCMlEjacLoO1P2yj66qHsBpOXOffjPWiO5t7WiIi4iWFEgkLyzYUcOW0TyiZN5oWZQfY7kohM28oyzbub+6piYiIlxRKJORVlZgfUTqX/lHfctyIZVzFw+wsspwoMS8iIkFPoURCWlWJ+f7WtTwSswiAJyvuJs84s7rE/KQlm3C6gqrvpIiIeKBQIiEtd2chcUU7eTHmVQDeqhzIItdV1e8bQIG9jNydhc00QxER8ZZCiYS0I0eP8FrMCyRYjvG1qztTKu/wOK6+UvQiItL8dEuwhC7D4LLvf0+SdS8HjFaMczxMRS1/pesqRS8iIsFBoURC1xd/JmnXUiqIZpxjAodoXWOIBUhOtJHRpU3Tz09ERHyi0zcSmravhE8mArClzxOsNbrX6H1T9Xzi4J5EWU99V0REgo1CiYSeQ1tg/i/BcMEFt3HekEfdSsxXSU60MXNkOlm9UppnniIi4hOdvpHQcqwQ5t4C5XbodCkMfhEsluoS87k7CzlYXEb7ePOUjY6QiIiEDoUSCR3OCpg/Cgp3QGInuOVdiI6rfjvKaqFv17bNOEEREWkMnb6R0GAY8K9fw4//wRl9Gp+m/4kvD1hVFE1EJIzoSImEhtWz4Ju/4cLCvcfuZ8VHx4GvSEm0MXFwT103IiISBnSkRILfpg8wlmUDMLXiVla4Lqp+a7+9TP1tRETChEKJBLddX2IsGoMFg3crBzDbeb3b2+pvIyISPhRKJHgdyoO5v8DiLOffznServwl1KhGov42IiLhQqFEglPxfpgzHMp+orB1b8ZXjMdJVJ0fUX8bEZHQplAiwafMDu8OB/tuaNOV7ZlvUEZcvR9TfxsRkdCmUCLBxVEK794M+7+H09rByEWkn3s2KYk2DyduTBYgRf1tRERCnkKJBI+KMph3O+z5CmyJMPJ9aNOFKKuFiYN7AjWvKFF/GxGR8KFQIsHBWQEL74Ydn0LMaXD7IkjpXf12Vq8U9bcREQlzKp4mzc/lhH/cD3lLISoObp0LnS6pMUz9bUREwptCiTQvlxMWj4cNC8EaDbe8A2ddXWOY02W4hZEbeqcqjIiIhBmFEmk+zkr45/3w/QKwRMHPX4fu/1dj2LINBUxasokC+4lbflVeXkQk/OiaEmkezgpYdLcZSKzRMOJNOG9YjWHLNhQwds5at0ACKi8vIhKOFEqk6VWWw/xRsOkDiIqFW+ZAz6E1hjldBpOWbMJT8XiVlxcRCT8KJdK0HKUw91bI+xdE2+AXc6HHII9Dc3cW1jhCcjKVlxcRCS+6pkSaTulheHcE5K+F6BZw2zw4q1+tw70tG6/y8iIi4UGhRJpG4U6Y83Mo3A4tWsNtCzze9nsyb8vGq7y8iEh40OkbCbz8dfDXa81Akngm3L283kACkNGljcrLi4hEEIUSCawNi+CNQVB6EM7oBaOXQ7vuXn1U5eVFRCKLQokEhssFK/9glo6vPA5nD4S7/gUJvtUVUXl5EZHIoWtKxP/K7PDBA/DDEvP55eMhcxJYoxq0OZWXFxGJDAol4l8F38L8O+HoTrMGyQ0vQp/bG73ZKKuFvl3bNn5+IiIStBRKxD8MA9a8AcuywVmOkdiJDZf/iR2WHrTffkRHNkREpF4KJdJ4pYfhwwnVp2sOpvTnjiO/JO8fx4H1gHrViIhI/XShqzTO5qXw6mVmILFGs/n833DpznvIK4pxG6ZeNSIiUh+FEmmY0sPwj/th3m1Qegja98Q5egV3bemL4aGyiHrViIhIfRRKxDcuF3zzN/jzRfDtXMACP3sY7s0ht6yTetWIiEiD6ZoS8d6+b+Cjx2Fvrvn8jPPhhuehUwYAB4sPe7UZ9aoRERFPFEqkfke2w4pnYNM/zeexp0P/JyDjPog68VdIvWpERKQxFEqkdkd3wecvwtq3wVUJWOCCX8A1T0FihxrDq3rV7LeX4emqEQtmJVb1qhEREU8USqSmA5vMMPL9QjCc5mvdroUBEyG5V60fq+pVM3bOWizgFkzUq0ZEROqjUCImZ4V5e+83b8KOnBOvn9Ufrvo1pF3h1WaqetVMWrLJ7aLXZNUpERGReiiURDLDgP3fm518v50LJQf+94YFzh0MVz4KqX183qx61YiISEMolIQ4p8vw7Ze/ywUF6yHvI9j4PhzZduK909pD+ijzp3XnRs1LvWpERMRXPoeS0tJSHnvsMVq3bk1JSQnTp08nLi6uxrjFixezYsUKysvLGT58OJmZmX6ZsJywbENBjdMkNcq5G4YZPPbkmqdltq+AY0dObCQqDrpfC+ePgB7XQZR7JVYREZGmYjEMw6fymqNGjWLYsGEMGzaMt99+m/Xr1/P888+7jdm8eTOjRo1i9erVGIbBxRdfzJIlS+jQoeYdG6cqKioiMTERu91OQkKCb6uJIMs2FDB2zlq3i0mjqSTNcoCzLfv41YWVdHPkwd6voewn9w/HxsNZV8O5Q6DHILDpf2cREWkcf/z+9imU5Ofn07VrV44ePYrNZuPQoUN07tyZAwcOEB8fXz3u/vvvp3379jzzzDMAPPDAA7Ru3ZopU6bU+x0KJR4YBjhK4fhROH4UZ2khU+atoMXx/aRYCkmxHKGz5SBplv3EWJw1Px9tM68NObMvnJ1pFjvTEREREfEjf/z+9un0TU5ODklJSdhsZvGrdu3aERsbS25uLgMGDKget3LlSn77299WP+/WrRuLFi3yuM3y8nLKy8urnxcVFfm0AK8dyoM1b5i/4KuOLzTqMY3fTlUeNFxQWQ6VZeZPxXH3x8d/AldF9VKigIkAHnJFqRHHNqMD24xU0i8bQJcLrzYrr0bHev+/lYiISDPwKZTs27ePNm3cC1/Fx8eTn59f5zhPY6pMnTqVSZMm+TKNhrHvgdWzAv89gRQVCy1aU2Q5ne9+slFgtKWANhQYbdlrtGObqwMFtKGqKshLHS6kixenzERERIKBT6HEYrFUHyWp4nA4iImJqXOcpzFVsrOzefTRR6ufFxUV0alTJ1+m5Z3WXeDKX1XNECwWHx/TiM/W89hiheg48zRLtA1iWrj/t0UraNEaYlqCxcLG7UcYOfurepescu4iIhJKfAolqamp2O12t9dKSkpITU2tc1xxcXGNMVXi4uI83r3jd227woCnA/89TUDl3EVEJBxZfRncr18/9u7di8PhAKg+JZORkeE2bsCAAWzZsqX6+bZt2+jfv39j5yr/U1XOHU46gPM/KucuIiKhyqdQkpqaSlZWFqtWrQJg+fLljBs3DpvNxowZM8jLywPMu28++eQTACorK8nNzWXMmDF+nnpkqyrnnpzofoomOdHGzJHpKucuIiIhx+c6JYcPH+bxxx8nLS2NwsJCpk2bRmxsLBdddBHZ2dkMHz4cgDfffJONGzficDgYNmyY10dKdEuwb3yu6CoiIhIATV6npCkolIiIiIQef/z+9un0jYiIiEigKJSIiIhIUFAoERERkaCgUCIiIiJBQaFEREREgoJCiYiIiAQFhRIREREJCgolIiIiEhQUSkRERCQo+NQluClUFZgtKipq5pmIiIiIt6p+bzemUHzQhZLi4mIAOnXq1MwzEREREV8VFxeTmJjYoM8GXe8bl8tFfn4+8fHxWCz+bSxXVFREp06d2LNnT1j31YmUdULkrFXrDD+RstZIWSdEzlprW6dhGBQXF5OamorV2rCrQ4LuSInVaqVjx44B/Y6EhISw/gtTJVLWCZGzVq0z/ETKWiNlnRA5a/W0zoYeIamiC11FREQkKCiUiIiISFCIqFASFxfHxIkTiYuLa+6pBFSkrBMiZ61aZ/iJlLVGyjohctYayHUG3YWuIiIiEpki6kiJiIiIBC+FEhEREQkKCiUiIiISFBRKREREJCgEXfG0higtLeWxxx6jdevWlJSUMH36dI9XBS9evJgVK1ZQXl7O8OHDyczM9Oq9YNHYddrtdu655x4+/vhjunbtysyZM7nsssuaehle8cefaZV58+Yxa9YscnJymmDmvvHXOo8dO8Zrr71G27Zt6d69e1D+ufpjrc8++ywulwur1cqxY8eYPHmy3ys/N5a36zx8+DB//OMfKS0t5ZVXXnF7L5z2R7WtMxz3R3X9mVYJh/1Rfets1P7ICAN33HGH8f777xuGYRhvvfWW8cgjj9QY88MPPxiXXHKJ4XK5DKfTafTp08fYu3dvve8Fk8auMzs723jvvfeMNWvWGIMHDzbat29vlJSUNOkavNXYtVbZt2+f0atXL+Pqq69uimn7zB/rPHLkiHHdddcZO3bsaLJ5N0Rj17pkyRLjnnvuqR571113GfPmzWuayfvAm3UahmF88803xpAhQ4w777zT7fVw2h8ZRu3rDLf9kWHUvtYq4bA/Moy619nY/VHIh5J9+/YZNpvNOH78uGEYhnHw4EGjRYsWRlFRkdu4++67z3jqqaeqn48bN8743e9+V+97wcIf61yxYkX163a73YiJiTFyc3ObYPa+8cdaq4wfP9545ZVXgnIn4K91Dhw40Pj000+bZM4N5Y+1/vGPfzRuueUWt/dee+21Jpi997xdZ5Wnnnqqxo49nPZHVTytM9z2R1U8rbVKOOyPqtS2zsbuj0L+mpKcnBySkpKw2WwAtGvXjtjYWHJzc93GrVy5ks6dO1c/79atG6tWrar3vWDhj3Vec8011a9X9SwIdJ+hhvDHWgFmz57NyJEjadmyZdNM3Ef+WOeHH37I1q1byc3N5brrriM7O5uKioqmW4SX/LHWIUOGsHjxYhYuXMjRo0c5fPgwd9xxR9MtwgverrOKp6Zl4bQ/quJpneG2P6pSWyO6cNkfVfG0Tn/sj0L+mpJ9+/bRpk0bt9fi4+PJz8+vc9zJY+p6L1j4Y50n27JlC/369SMlJSUwE24Ef6x127ZtFBUVkZGRwaZNmwI/6QbwxzrfffddLr30UiZMmMDo0aPp06cPhmEwbdq0wC/AB/5Ya48ePXjvvfe47bbbuPrqq1m0aBEtWrQI/OR94O06fdlGKO+PvBUO+6O6hNP+qC7+2B+F/JESi8VSneyqOBwOYmJi6hx38pi63gsW/ljnyV555RWee+65wEy2kRq7VqfTyauvvsqECROaYroN5o8/040bN3LFFVcQGxtL27ZtGTNmDG+//XbgJ+8jf/39LSsrY/HixWzfvp1bb72VysrKwE7cR96u05dthPL+yFvhsD+qTbjtj+rij/1RyIeS1NRU7Ha722slJSWkpqbWOa64uLh6TF3vBQt/rLPKxx9/zDXXXEOXLl0CN+FGaOxav/jiC2bNmkXbtm1p1aoV48aN47///S+tWrVqiul7zR9/ppWVlTidzur3evfuTWFhYQBn3TD+WGtOTg5r1qwhMzOTnJwc1q1bx/PPPx/4yfvA23X6so1Q3h95I1z2R7UJt/1RXfyxPwr5UNKvXz/27t2Lw+EAqD7UlJGR4TZuwIABbNmypfr5tm3b6N+/f73vBQt/rBPMJLt7926GDh3aBLNumMau9ZJLLmHTpk2sX7+e9evX88wzz3DxxRezfv36JluDN/zxZ9q7d2+2bt1a/V50dDTnnHNOoKfuM3+sdeHChXTr1g0wd6BTpkzhs88+a4rpe83bddYlnPZH9Qmn/VFtwm1/VBd/7I9CPpSkpqaSlZVVfSHY8uXLGTduHDabjRkzZpCXlwfA/fffzyeffAKYaS43N5cxY8bU+16w8Mc6t2/fzuzZsxk4cCA//vgj3377LTNnzmyeBdWhsWu12WykpaVV/1RdvJWWltZcS/LIH3+mDz/8MEuXLqW8vBww/1U2fvz4ZlhN3fyx1gsvvJB169ZVb9Nisfj8SzDQvF1nFcO8A9LttXDaH1XxtM5w2x9VOXWt4bY/quLpz9Qv+6MG37cTRA4dOmSMHj3amDx5svHII48Y5eXlhmEYRnp6urFgwYLqcW+88Ybxq1/9yhg/fryxcuVKt23U9V6waMw68/PzjY4dOxqA208w1nkwDP/8mVZ58803g/IWPMPwzzrfeecdY8yYMcaMGTOMZ555pknn74vGrtXpdBpPPfWUMXXqVGP27NnG008/bZSVlTX5Ourj7TrXrFljXHnllcY555xjrFq1ym0b4bQ/8rTOcN0f1fVnWiUc9kd1rbOx+yOLYZwSdURERESaQcifvhEREZHwoFAiIiIiQUGhRERERIKCQomIiIgEBYUSERERCQoKJSIiIhIUFEpEREQkKCiUiIiISFBQKBEREZGgoFAiIiIiQUGhRERERIKCQomIiIgEhf8P1c/P+P+39nMAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "x = np.linspace(0,0.15,500)\n",
    "plt.plot(w,pcw,'o')\n",
    "plt.plot(x,lognorm.cdf(x,paramsw[0],scale=paramsw[2]))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "ff4b90f4-10e9-4e59-ad7c-6af7acdb5d1d",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "be9d268c-9c20-4341-a820-b4326fcd4bbb",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "1ce6e482-e9c8-4202-a8c9-d94496d024d6",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.11.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
